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PREFACE 


The  purpose  of  this  study  was  to  develop  an  alternative  systeti  to 
the  present  plan  for  evacuating  and  distributing  casualties  resulting 
fron  a  conventional  European  conflict.  The  proposed  systen  requires 
Modified  Boeing-767  aircraft,  as  part  of  the  Civil  Reserve  Air  Fleet 
( CRAF ),  to  perforn  the  intertheater  evacuation  Mission.  The  C-9 
Nightingale  aircraft  uould  distribute  casualties  fron  a  single  B-767/C-9 
interface  location.  The  plan  relieves  the  C-141  Starlifter  of  all 
evacuation  conMitnents  so  that  it  May  be  dedicated  to  the  resupply 
effort. 

I  wish  to  thank  Many  people  who  have  contributed  so  Much  to  this 
research.  These  include  Br.  Jeffrey  Kennington  ( Departneni  of 
Operations  Research,  Southern  Methodist  University,  Dallas,  Texas),  Mrs. 
Georgia  Thonas  ( Departnent  of  Health  and  HuMan  Services,  Rockville, 
Maryland)  and  TSgt.  Linnes  Chester  (Medical  Readiness  Division,  Fort  SaM 
Houston,  Texas).  I  an  very  grateful  to  Capt.  Robert  ChMielewski  and 
Maj.  Robert  Murray  at  MAC'S  Operations  Research  Division  and  Surgeon 
General's  Office,  respectively,  for  their  tiMe,  research  data  and 
recoMMendations  for  this  study.  Their  insight  greatly  increased  My 
understanding  of  the  problens  of  this  real-world  possibility. 

Special  thanks  are  extended  to  Capt.  Joseph  Litko,  mv  thesis 
advisor,  who  was  always  available  to  guide  Me  through  those  difficult 
tines.  His  wisdoM,  encouragenent  and  recoMMendations  Made  this  thesis 
possible.  Thank  you,  Joe. 


Uillian  B.  Ewing,  Jr. 
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^^The  objective  of  this  FTrsarch  was  to  develop  and  evaluate  a 
casualty  evacuation  and  distribution  system  using  B-?6?  and  C-9  aircraft. 
For  a  European  conventional  conflict,  an  average  casualty  rate  between 
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The  performance  of  the  system  was  measured  by  the  average  time  each 
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the  patient  arrived  at  the  final  destination  airport.  The  factors  in 
the  model  affecting  the  mean  time  in  the  system  (TIS)  include  pre- 
departure  ground  time,  flying  time,  number  and  capacity  of  aircraft  and 
casualty  rate. 

Response  surface  equations  were  generated  from  the  experimental 
results  for  selected  combinations  of  factor  levels.  The  prediction 
equations  provide  an  accurate  measure  of  the  performance  of  the  system, 
while  saving  the  time  and  expense  of  conducting  simulation  experiments. 
The  equations  can  be  used  to  determine  either  the  required  number  of 
aircraft  or  the  necessary  aircraft  capacity  given  a  specified  criterion 
value  of  mean  TIS  and  expected  casualty  rate.  Kc- u * 
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ABSTRACT 

The  objective  of  this  research  was  to  develop  and  avaluata  a 
casualty  avacuation  and  distribution  systan  using  B -767  and  C-9 
aircraft.  For  a  Europaan  convantional  conflict,  an  avaraga  casualty  rata 
batwaan  1600  and  1900  par  day  was  considarad  for  a  60- day  pariod. 
Casualtias  wara  distributad  anong  all  potantial  nenbers  of  tha  National 
Disastar  Madical  Systan.  Tha  systan  was  nodal lad  using  SLAM  sinulation 
and  FORTRAN  conputer  coda. 

Tha  perfornance  of  tha  systan  was  Measured  by  tha  average  tine  each 
patient  spent  in  tha  evacuation  systan,  beginning  with  the  tine  tha 
patient  was  released  (nedically  cleared)  to  fly.  Tha  tine  ended  when 
tha  patient  arrived  at  tha  final  destination  airport.  Tha  factors  in 
tha  nodal  affecting  tha  naan  tine  in  tha  systan  ( TIS )  include  pre¬ 
departure  ground  tine,  flying  tine,  nunbar  and  capacity  of  aircraft  and 
casualty  rata. 

Response  surface  aquations  wara  generated  fron  tha  experinental 
results  for  selected  conbinations  of  factor  levels.  Tha  prediction 
aquations  provide  an  accurate  neasure  of  tha  parfornanca  of  the  systan, 
while  saving  tha  tine  and  expanse  of  conducting  sinulation  exparinents, 
Tha  aquations  can  be  used  to  deternine  either  the  required  nunber  of 
aircraft  or  tha  necessary  aircraft  capacity  given  a  specified  criterion 
value  of  naan  TIS  and  expected  casualty  rata. 
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CASUALTY  EVACUATION  AND  DISTRIBUTION  USING  B -767  AND  C-9  AIRCRAFT 


I.  CONCEPTUALIZATION 


iDlEQdUCiiQD 

A  great  challenge  facing  the  Military  Airlift  CoMMar.d  (HAC)  is  to 
satisfy  the  airlift  requirenents  of  the  Dcpartnent  of  Defense  (DoD)  in 
the  event  of  a  Major  conventional  war  in  Europe  and  other  regional 
conflicts  at  the  sane  tine.  One  inportant  airlift  requirenent  that  this 
research  effort  investigates  is  the  European  evacuation  and  CONUS 
distribution  of  wartine  casualties. 

In  a  Congressionally  Mandated  Mobility  Study  (CMMS)  in  1981,  the 
DoD  deternined  that  an  airlift  capability  of  66  Million-ton-niles  per 
day  (MTM/D)  was  necessary  to  support  all  conbat  forces  (10i53).  Today, 
MAC  has  about  a  35  MTM/D  capability;  but  in  referring  to  the  66  MTM/D 
goal,  General  Thonas  Ryan,  CINCMAC,  stated,  "...  we  will  attain  it 
before  the  end  of  the  century"  (10>53).  Recent  airlift  additions  and 
Modifications  are  closing  the  gap. 

EobdDCeaeota.  MAC  has  Made  considerable  progress  toward  a  66  MTM/D 
capability.  As  reported  in  the  Novenber  issue  of  Air  Force  Magazine 
<10a33-89),  airlift  enhanceMents  include  all  MAC  forces.  Fi.st,  the  C~ 
141  "Starlifter"  fl»,et  has  been  Modified  for  aerial  refueling  and  its 
cargo  conpartnent  has  been  stretched  thus  increasing  its  bulk  cargo 
capacity  by  30X.  The  C-5A  "Galaxy"  is  being  Modified  with  stronger  wings 
and  engines.  Also,  the  first  of  MAC'S  newest  airlift  addition,  the  C-5B, 
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will  be  delivered  by  DeceMber  1985  and  the  fiftieth  thould  be  delivered 
in  1989.  Additionally,  full  scale  development  of  the  C-17  was  approved 
in  the  summer  of  1985.  If  fully  funded,  210  C-17  aircraft  will  be 
produced.  The  C-130  "Hercules"  aircraft  outer  wing  boxes  are  being 
replaced,  extending  its  flying  life  up  to  40,000  hours.  Also,  a  joint 
SAC-HAC  procurement  of  40  KC-10  aircraft  will  enhance  airlift  as  an  in¬ 
flight  refueler  and  a  cargo  plane.  Finally,  some  aircraft  of  the  Civil 
Reserve  Air  Fleet  (CRAF)  are  being  Modified  to  carry  Military  cargo(10). 

dtE&QCdiGll  CBflE-  Airlift  augMentation  by  U.S.  owned  coMMercial 
aircraft  will  be  iMpleHented  in  stages  depending  on  the  severity  of  any 
national  emergency.  The  CRAF  will  transport  soldiers  and  supplies  to 
the  war  theater  and  May  also  be  tasked  to  evacuate  Military  and  state 
departnent  dependents.  MAC  has  considered  Modifying  the  Boeing  aircraft 
( B-767 )  during  its  production  phase  to  Make  it  aeroMedically  capable 
(27).  The  B-767  would  be  used  prinarily  for  intertheater  evacuation, 
while  the  C-9  "Nightingale"  would  be  the  priMary  intratheater  evacuation 
aircraft.  Smaller  CRAF  aircraft  could  also  be  Modified  to  someday 
replace  the  C-9  if  they  are  required  for  the  tactical  evacuation  of 
casualties.  These  Modified  CRAF  would  be  incorporated  into  the  CRAF 
plan  and  used  when  necessary. 

Accaacdical  Beauirfifler/i .  The  Medical  Readiness  Division  of  the 
Office  of  the  Surgeon  General  anticipates  the  saturation  of  U.S. 

Military  hospitals  in  Europe  in  the  early  stages  of  the  conflict.  As  a 
result,  the  DoD  plans  to  evacuate  all  casualties  that  cannot  be  returned 
to  duty  within  15  days  (17). 


Under  current  plant,  the  C-141  will  be  used  for  tha  strategic 
evacuation  of  thaia  catualiiat.  Thit  will  ba  acconplithed  during  tha 
raturn  (backhaul)  portion  of  iti  ratupply  Minion.  Tha  nadical  Mittion 
could  inpact  tha  ratupply  Mittion  but  iti  iMpact  it  difficult  to 
attiMata.  An  exanination  of  tha  extra  tiMa  required  for  a  C-141  to 
execute  a  coMbinad  ratupply  and  avacuation  Mittion  May  ba  inatructiva. 

Btiupply  list*  Tha  avacuation  iMpact  on  tha  ratupply  Mittion  will 
ba  axaHinad  in  terns  of  tha  differanca  in  turn-around-tines  for  a  C-141 
Mittion  with  and  without  tha  avacuation  requirenent.  The  planned 
MaxiMUM  ground  tine  for  cargo  onload  or  offload  it  two  hourt  and  fifteen 
Hinutet  (23).  Concurrent  refueling  operation!  and  optiwal  crew  ataging 
are  attuned  for  the  contingency.  Therefore,  alapted  turn-around-tine 
thould  not  exceed  four  and  one-half  hourt  plut  tha  enroute  flying  tinea 
for  the  C-141  retupply  nittion. 

With  the  nadical  evacuation  taeking,  retupply  turn-around-tine  will 
include  the  following  additional  activitiea  and  eatiMated  average  tinea 
to  conplete  each  activity  '&)> 

1.  Repoiitioning  to  caaualty  onload  location  -  45  Minutes 

2.  Reconfiguring  aircraft  for  casualties  -  45  Minutes 

3.  Uploading  casualties  at  onload  location  -  90  Minutes 

4.  Offloading  casualties  at  offload  location  -  90  Minutes 

5.  Reconfiguring  aircraft  for  cargo  -  30  Minutes 

These  activities  and  estinates  are  based  on  the  67  litter  patient 
capacity  of  the  C-141  (23).  Optinal  crew  staging  is  still  assuned.  The 
elapsed  turn-around-tine  for  the  Medical  Mission  is  five  hours  longer 
than  that  for  the  resupply  Mission  alone.  The  possibility  of  a  Medical 
diversion  was  not  considered  but  a  nininun  additional  one  hour  delay 


would  be  incurred. 


I 
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In  addition  to  tha  resupply  inpact,  othor  concerns  should  bo 
Mentioned.  Hospital  bods  or  C-9  aircraft  support  mist  bo  availablo  at 
tha  offload  location.  This  is  required  bocauso  tha  offload  location  has 
already  boon  deternined  by  the  rosupply  requirenents  and  cannot  ba 
adjusted  for  patient  convenience.  Also,  the  nodical  oquipnont 
(stanchions,  litters,  etc.)  returning  to  Europe  will  reduce  tha  anount 
of  prinary  resupply  cargo  that  can  ba  carried  on  each  aircraft. 

The  tine  factor  alone  is  insufficient  to  deternin#  tha  inpact  each 
nodical  Mission  has  on  the  resupply  nission.  The  C-9  scheduling  and  the 
reduced  anount  of  prinary  war  supplies  carried  on  C-141  aircraft  are 
factors  that  nust  also  be  considered.  The  inpact  will  be  reduced  or 

i 

elininated  under  the  proposal  studied  in  this  research. 

I 

Ecapasal 

This  research  effort  will  exanine  an  evacuation  and  distribution 
systan  using  B -767  and  C-9  aircraft.  The  b-767  aircraft  will  be  used 
prinarily  for  the  strategic  evacuation  Mission,  while  the  C-9  aircraft 
will  be  used  for  the  i.itratheater  evacuation  nission  in  the  CONUSf 
however,  any  Modified  aircraft  could  execute  the  stateside  nission.  If 
the  C-9  is  deployed  to  Europe  for  tactical  evacuation  of  casualties,  a 
replacenent  aircraft  ( probably  CRAF )  would  be  required. 

Bcscaccb  Question 

As  an  alternative  to  the  present  evacuation  plan,  what  is  the 
structure  of  an  evacuation  and  distribution  network  using  Civil  Reserve 

/ 

Air  Fleet  <  CRAF )  and  C-9  aircraft  resources  and  how  can  its  perfornance 
be  evaluated? 
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BackQEfiuod 

The  existing  plan  and  any  alternatives  addressing  the  avacuation 
and  distribution  of  casualties  it  an  intansiva  routing  and  scheduling 
problen.  It  uas  originally  investigated  by  Lt.  Col.  McLain  at  HAC's 
Operations  Research  Division  as  a  nulti-connodity  network  flow  problen 
(20);  his  Casualty  Evacuation  Model  (CEM)  included  the  followingi 

1.  9  European  onload  locations. 

2.  73  CONUS  offload  locations. 

3.  95  enroute  cargo  terminals. 

4.  11  patient  categories. 

5.  60  one-day  tine  periods. 

A  review  of  the  literature  pertaining  to  this  kind  of  problen  will 
begin  with  the  classical  approach  to  solving  network  related  problens. 

A  very  recent  study  that  exanined  sone  variations  of  the  original 
problen  will  also  be  reviewed.  Finally,  the  latest  study  on  solving 
large  network  problens  will  be  exanined. 

Blit  Belated  Studies.  The  casualty  evacuation  and  distribution 
problen,  notwithstanding  its  size,  is  a  classical  network  problen  if  the 
patient  categories  and  the  tine  periods  are  disregarded.  In  this  case, 
one  attenpts  to  find  the  shortest  path  o»  naxinun  flow  solutions.  If 
the  tine  periods  are  included,  the  problen  becones  a  nulti-period 
transportation  or  transshipment  problen  with  the  objective  to  nininize 
shipping  costs.  Each  approach  will  be  exanined. 

dtlUQCk  pcableo.  The  casualty  evacuation  problen  has  well- 
defined  sources,  sinks  and  a  connected  systen  of  arcs  on  which  a 
connodity  flows  (31).  Docunentation  for  solving  the  network  for  naxinun 
flow  or  shortest  path  includes  that  provided  by  Hillier  and  Liebernan 
(15:241)  or  Uu  and  Coppins  (36). 


To  MaxiHize  flow,  the  three-step  algorithm  by  Hilliar  and  Liebernan 
night  ba  used  but  tha  routa  structure  Must  alraady  ba  established,  in 
which  casa  tha  problen  can  ba  raducad  to  assigning  aircraft  to  routas 
<S>.  Unfortunately,  tha  CEM  is  not  static}  tha  supply  and  danand  at 
tha  sourcas  and  sinks  changa  daily,  and  tha  routa  structura  has  not  bean 
astablishad.  As  a  result,  other  techniques  will  ha  reviewed. 

Tha  four-step  algorithm  presented  by  Wu  and  Coppins  can  be  used  to 
find  tha  shortest  paths  in  tha  network.  A  research  effort  by  Cuanthar 
involved  tha  distribution  of  Hail  to  six  cities  located  within  50  Miles 
of  each  other  (13).  In  addition  to  finding  tha  shortest  path,  ha  showed 
that  exhaustive  emmeration  (complete  identification)  of  all  possible 
routas  would  result  in  an  optiMal  solution.  Unfortunataly,  with  larger 
networks  ( ha  triad  a  10-node  systen ),  tha  exhaustive  enuneration 
technique  becane  prohibitively  expansive.  Tha  underlying  conplication 
in  tha  CEM  is  tha  size  of  tha  problen. 

Another  shortest  path  approach  investigated  by  Guenther  had  also 
bean  used  by  Gaskall  (12).  By  using  various  refining  heuristics  (a 
rule-of-thunb  for  solving  sone  aspect  of  a  probleM  (7*623)),  Guenther 
avoided  tha  conplete  enuneration  of  feasible  solutions  by  early 
elinination  of  partial  solutions  known  to  produce  sub-optinal  final 
solutions.  If  tha  CEH  is  sub-divided  into  a  systen  of  snaller  networks, 
this  approach  could  be  used  to  select  a  shortest  path  routing  structure 
for  each  subsysten.  Unfortunately,  a  nore  difficult  problen  arises  - 
what  criteria  should  be  used  for  sub-dividing  the  network  and  then,  in 
what  sense  is  optinality  defined  in  the  larger  systen? 
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IcaDlBfidaiiQO  EcflblttJ.  The  CEM  is  also  sinilar  to  a  special 
category  of  network  problens  called  transportation  problems.  Instead  of 
Haxinizing  flow  or  nininizing  distance,  the  transportation  probleti  has 
associated  shipping  costs  for  each  route  between  source  and  sink.  The 
objective  is  to  choose  the  conbination  of  routes  that  nininizes  total 
shipping  or  distribution  costs  (32H87).  Hillier  and  Liebernan  also 
provide  Methods  for  solving  transportation  problens  (15H34). 

In  separate  research  on  the  Non-conbatant  Evacuation  Operation 
(NEO),  Gullett  (14)  and  Moncure  (24)  observed  that  the  size  of  the 
problen  prohibited  the  application  of  any  of  the  Methods  presented  by 
Hillier  and  Liebernan.  The  NEO  deals  with  the  evacuation  of  Military 
and  State  Departnent  dependents  fron  Gernany  in  the  event  of  a  conflict. 
A  one-day  snapshot  of  the  CEH  without  regard  for  the  patient  categories 
would  require  a  square  Matrix  with  177  rows,  or  less  than  100  rows  if 
the  nuMber  of  interMediate  transfer  locations  could  be  reduced  to  less 
than  a  dozen  (15U50). 

The  solutions  would  be  one-day  estiMates  of  a  particular  static 
scenario  rather  than  the  actual  dynanic  evacuation  process.  This 
liMitation,  together  with  the  conplexity  of  solving  a  nulti-coMModity 
flow  problen,  substantially  reduces  the  applicability  of  classical 
network  approaches  to  solving  the  casualty  evacuation  problen. 

Becenl  Study*  The  groundwork  for  the  Casualty  Evacuation  Model 
(CEM)  was  acconplished  by  Lt.  Col.  McLain,  Deputy  Chief  of  the 
Operations  Research  Division,  MAC  Headquarters.  In  March  1985,  Alfano 
and  O'Neill  conpleted  their  study  of  a  Modified  CEM  < 1 ).  They  assunod 
that  all  casualties  would  arrive  at  a  centrally  located  eastcoast 


Military  basa.  Fro*  that  location,  CRAF  aircraft  would  transport 
pationts  to  two  othar  distribution  cantors  in  tha  CONUS.  C-9  aircraft 
would  distributa  up  to  40  patiants  (35H56)  to  tha  various  hospitals 
froM  ona  of  tha  thraa  distribution  cantors.  An  avaraga  of  1000 
casualtias  par  day  was  usad  in  thair  Modal. 

Tha  hub-and-spoka  Modal  davalopad  by  thaw  was  tastad  in  thraa 
phasas.  Tha  first  allowad  an  unliwitad  nuMbar  of  C-9  aircraft.  It 
raquirad  up  to  17  C-9s  and  tha  avaraga  patiant  spant  3.3  hours  in  tha 
natwork.  Tha  sacond  phasa  liHitad  tha  nuMbar  of  C-9s  in  aach 
distribution  araa  to  achiava  tha  sana  avaraga  patiant  tiMa.  Tha  Modal 
raquirad  at  laast  7  aircraft  par  distribution  cantor.  Tha  last  phasa 
allowad  patiants  to  quaua  at  tha  distribution  cantors  in  ordar  to  fill 
aach  C-9.  This  varsion  raquirad  10  C-9  aircraft  (U78). 

Tha  final  conclusion  frow  Alfano  and  O'Naill's  study  was  that  an 
al tomato  distribution  systaw  should  ba  davalopad  to  raduca  tha  C-9 
raquiranant.  Tha  raason  for  this  -  tha  Air  Forca  has  only  11  CONUS 
basad  C-9  aircraft  (1>75). 

Cucctui  Status.  OptiMality  in  tarns  of  aithar  tha  shortast  path  or 
waxiMUM  flow  for  a  natwork  of  this  siza  has  navar  boon  achiavad  dua  to 
tha  problaw's  cowbinatorial  natura  and  thus,  "beyond  tha  scopa  of 
axisting  coxputar  codas"  (20).  Howavar,  continuad  intarast  in  this 
problaM  has  prowptad  tha  Oparations  Rasaarch  Departnant  of  Southarn 
Methodist  Univarsity  (SNU)  to  subwit  a  proposal  to  tha  Air  Forca  Offica 
of  Sciantific  Rasaarch  for  a  grant  to  davalop  and  avaluata  a  casualty 
avacuation  Modal  for  a  European  conflict  ( 17 ).  Tha  grant  was  Made  and 
tha  study  should  bagin  in  OacaMbar  1985?  howavar,  praliMinary  results 
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of  a  snaller  version  of  tho  original  problen  war*  presented  by  Dr.  Allan 

i 

of  tha  SHU  staff  (2). 

Dr.  Allan's  constraint  Matrix  of  McLain's  original  varsion  bad  More 
than  137,000  rows.  To  solva  this  problen,  Dr.  Allan  plannad  to  davalop 
a  spacializad  procadura  to  solva  a  scalad  down  varsion  of  tha  original 
problen.  Tina  periods,  patiant  catagorias  and  facility  locations  could 
ba  aggragatad  to  raduca  tha  siza  of  tha  problaM. 

A  procadura  was  davalopad  to  solva  a  tast  problaM  of  about  9000 
rows.  It  involvad  ganarating  uppar  bounds  on  tha  objactiva  function  by 
solving  tha  problaM  using  a  rasourca-diractiva  dacoMposition  technique 

t 

(151704),  and  than  ganarating  lowar  bounds  on  tha  objactiva  function  by 
partially  solving  a  dual  of  tha  problaM  (2>8).  Tha  itarative  procasss 
stops  whan  tha  diffaranca  batwaan  both  bounds  ranches  a  cartain 

tolaranca.  Tha  conputer  prograM  "EVAC"  davalopad  by  Allan  raquirad  lass 

i 

MeMory  and  about  65X  lass  tiMe  than  tha  prograM  "MCNF"  davalopad  by 
Kannington  to  solva  a  sinilar  problaM  (2i48>.  Tha  diffarancas  wara  due 
to  Kannington' s  prograM  saaking  an  axact  optinun.  "EVAC"  should  ba  a 
good  Modal  for  tha  scalad  down  casualty  evacuation  problaM. 

Tha  SHU  Operations  Research  staff  have  also  proposed  to  investigate 
tha  applicability  of  using  tha  projective  algorithM  of  Dr.  Narandra 
Karnarkar  of  AT&T  Ball  Labs  to  selva  tha  full  nodal  (17).  A  tast  case 
using  an  SHU  varsion  of  Karnarkar's  algorithn  to  solva  assignnent 
probleMS  was  not  encouraging?  however,  tha  poor  result  nay  be  due,  in 
part,  to  an  inaccurate  version  or  application  of  tha  algorithn  (3). 


BftVitW  CfiOCllilifiD*  Many  probl an*  of  this  typo  have  baan 
solved  by  various  nethods;  however,  nona  of  tha  problens  approachad  tha 
siza  of  tha  CEH.  Qptinality,  though  desirable,  is  not  tha  goal  of  this 
study.  A  faasibla  solution  will  ba  found  and  avaluatad  to  provida  a 
basis  for  conparison  and  futhar  study  in  this  araa.  Tha  nethodology 
usad  by  Gullet,  Moncure  and  Alfano  will  also  ba  usad  in  this  rasaarch 
affort  -  sinulation  -  to  ba  discussed  in  Chaptar  II. 

Iaciical  Btepacatifio 

According  to  tha  DoD  Thaatar  Conbat  Hadical  Systen  (Figura  1 >, 
casualtias  will  ba  procassad  through  four  achalons  of  nodical  cara 
bafora  raaching  thair  recovary  hospital. 
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Figure  1.  DoD 

Thaatar  Conbat 

Medical  Systen 

Tha  first  echelon  is  self-aid  or  buddy-care  in  the  field.  At  the 
second  echelon,  wounds  and  tha  general  conditions  of  each  patient  will 
be  evaluated  to  deternine  the  priority  for  further  treatnant  and 
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evacuation.  The  third  echelon  it  the  first  Medical  facility  staffed  and 
equipped  to  provide  specialty  care.  General  surgery  and  further 
evacuation  determinations  will  be  performed  here.  Comprehensive  medical 
care  will  be  provided  at  the  fourth  echelon  in  theater  if  the  patient 
can  be  rehabilitated  within  15  days.  If  not,  evacuation  to  the  CONUS 
will  be  coordinated  through  the  Armed  Services  Medical  Regulating  Office 
(ASMR0>  and  the  Joint  Medical  Regulating  Office  (8>. 

NtiiiQDal  Diiailcc  Utdical  Sysieu 

In  December  1781,  President  Reagan  established  the  Emergency 
Mobilization  Preparedness  Board  to  develop  national  policy  and  programs 
to  improve  emergency  preparedness.  A  National  Disaster  Medical  System 
(NGMS)  was  developed  in  response  to  the  President's  mandate.  The  NDMS 
task  force  includes  the  Department  of  Defense,  Department  of  Health  and 
Human  Services,  the  Federal  Emergency  Management  Agency  and  the  Veterans 
Administration  (29). 

The  NDMS  is  designed  to  care  for  the  casualties  of  any  incident 
that  exceeds  the  medical  care  capability  of  an  affected  region.  Two 
often  quoted  situations  for  which  NDMS  was  designed  include  a  major 
California  earthquake  and  a  conventional  Europwan  war.  NDMS  will  inclde 
all  major  metropolitan  areas  of  the  nation  (26).  It  incorporates  the 
Civilian-Military  Contingency  Hospital  System  (CMCHS)  that  already 
includes  770  hospitals  in  48  metropolitan  areas  (29).  The  additional 
proposed  metropolitan  areas  are  shown  in  Figure  2. 

H6C  UiaSifiD 

The  Military  Airlift  Command  will  provide  the  airlift  resources  to 
connect  the  operations  of  the  DoD  Theater  Combat  Medical  System  and  the 
National  Disaster  Medical  System. 
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NATIONAL  DISASTER  MEDICAL  SYSTEM 
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Btxtaccb  Liaiii 

This  research  effort  will  examine  on*  evacuation  and  distribution 
concept  starting  with  th*  patient's  r*l*as*  from  th*  European  fourth 
•ch*lon  of  car*  to  his  arrival  at  on*  of  94  U.S.  Metropolitan  area 
airports  (9).  Th*  operation  of  th*  evacuation  and  distribution  systeN 
will  extend  through  th*  first  60  days  of  th*  conflict.  Th*  estimated 
number  of  casualties  will  vary  between  1000  and  2500  per  day.  Other 
parameters  of  th*  Model  will  b*  discussed  later. 

Htfiiucei  q£  EfftcUytneis 

In  an  operational  environment  involving  expensiv*  and  linit*d 
resources,  reasonable  Measures  of  Mission  effectiveness  include  aircraft 
and  crew  utilization,  average  cargo  capacity  utilization,  and  the 
average  tine  required  to  transport  a  given  tonnage  to  a  specific 
location  ( closure  tine ).  Fron  a  Medical  perspective,  these  Measures 
would  be  superceded  by  th*  patient's  average  tine  in  the  systen, 
expected  environmental  conditions  and  th*  maximum  tin*  any  patient  Might 
spend  in  the  evacuation  and  distribution  system. 

To  varying  degrees,  all  these  Measures  will  be  considered  in  the 
proposed  Model  of  the  systen.  Th*  primary  Measure  of  effectiveness  will 
be  th*  average  tin*  a  patient  spends  in  th*  system  (response  variable). 
This  time  will  begin  when  a  full  load  of  patients  is  released  for 
evacuation  and  will  end  when  those  patients  arrive  at  their  final 
destination  airport. 

It  is  important  to  not*  that  a  large  average  time  in  the  system 
does  not  mean  that  a  patient  spends  that  much  time  travelling  in  an 
aircraft  or  bus.  The  response  variable  (time  in  the  systeM)  includes 
that  time  spent  in  the  hospital  waiting  for  transportation. 
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Secondary  measures  of  effectiveness  includa  rasourca  utilization 
rates,  maximum  time  in  tha  system,  and  tha  nunbar  of  patients  awaiting 
evacuation.  Alto,  tinea  all  CRAF  aircraft  will  depart  with  a  full  load 
of  casualties,  maximum  cargo  utilization  rata  will  be  achieved. 

Btaeiccb  Objective 

A  routing  network  between  a  casualty  onload  location,  intermediate 
locations,  and  an  aggregation  of  75  final  offload  locations  will  be 
designed  and  evaluated  under  various  conditions.  These  conditions 
include  relatively  high  and  low  values  of  the  Major  factors  affecting 
the  prinary  Measure  of  effectiveness.  The  factors  include  daily 
casualty  rate,  CRAF  capacity,  ground  tiMes,  flying  tine,  nunber  of  CRAF 
aircraft  and  nuMber  of  C-9  aircraft. 

Sub-Objectives.  Prinary  sub-objectives  include  detemining  how 
Many  inter Mediate  locations  will  be  required  and  where  they  should  be 
located.  Aircraft  patient  capacity  and  average  daily  casualty  rate  will 
be  determined  in  order  to  estimate  the  average  nunber  of  loads  departing 
Europe.  From  this,  an  estimate  of  the  required  number  of  CRAF  aircraft 
can  be  made.  With  the  simulation  model  of  the  network,  aircraft 
utilization  rates  will  be  determined. 

Subuecy 

Airlift  plays  a  vital  role  in  the  conduct  of  all  contingency 
operations.  Potential  airlift  shortages  during  a  conventional  European 
conflict  could  have  disastrous  consequences.  The  recent  force 
Modifications,  planned  additions  and  contingency  augmentation  plans  may 
not  be  sufficient  to  meet  all  the  DoD  wartime  demands,  including 
casualty  evacuation. 
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The  usual  techniques  for  solving  network  problems  are  not 
appropriate  for  the  CEH  because  of  its  large  size.  In  a  review  of 
related  studies,  this  limitation  was  recognized.  Simulation  techniques 
used  by  Alfano  and  O'Neill  seen  to  be  the  Host  appropriate  Method  for 
addressing  this  probleH  and  will  be  discussed  in  Chapter  II. 

This  study  will  examine  an  alternative  evacuation  and  distribution 
plan  based  solely  on  CRAF  and  C-9  resources.  It  will  begin  with  the 
patient's  release  in  the  European  fourth  echelon  of  care  and  end  with 
his  arrival  at  the  airport  nearest  a  NDHS  participating  hospital.  The 
primary  Measure  of  effectiveness  will  be  the  patient's  time  in  the 
system,  controlled  by  the  six  major  factors  included  in  the  network. 

This  performance  measure  will  be  used  to  evaluate  the  network  at 
different  factor  levels.  The  study  will  attempt  to  develop  an 
evacuation  and  distribution  network  that  results  in  acceptable  measures 
of  performance  including  the  average  and  maximum  time  in  the  system,  and 
resource  utilization. 
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XI.  MODEL  FORMULATION 


InicfidtedifiD 

The  casualty  evacuation  and  distribution  problen  is  a  conplex 
interaction  of  radical  requirenents,  aircraft  and  aircrew  scheduling  and 
support  operations.  It  is  a  dynamic  problen  that  has  optinal  or  near- 
optinal  solutions  for  each  situation  or  significant  tine  period  in 
question.  The  focus  of  this  research  effort,  however,  is  on  designing  a 
sinulation  Model  of  an  evacuation  and  distribution  network  and 
evaluating  the  Model's  per for nance. 

Simulation 

The  techniques  of  siMulation  are  well  suited  for  this  type  of 
probleH.  Banks  and  Carson  define  siMulation  as  "the  iHitation  of  the 
operation  of  a  real-world  process  or  systeri  over  tiwe"  (4«2). 

Sinulation  is  a  Method  a  research  analyst  can  use  to  evaluate,  predict 
and  soMetines,  optinize  the  operation  of  a  systeM  that  can  not  be 
fornulated  analytically  <  33120).  With  the  flexibility  provided  by 
sinulation,  one  can  experinent  with  the  systeM  to  deterHine  how  it 
reacts  to  changes  before  coMMitting  resources. 

For  these  reasons  and  others,  Alfano  and  O'Neill  enployed 
sinulation  techniques  to  evaluate  their  proposed  evacuation  and 
distribution  systeM  (U20).  They  sinulated  the  systeM  for  60  days 

beginning  with  the  arrival  of  1000  patients  per  day  at  an  eastcoast 

/ 

location  via  C-141  aircraft. 

Using  CRAF  and  C-?  aircraft  in  a  hub-and-spoke  systeM,  the  Model 
sinulated  the  distribution  of  casualties  anong  73  locations.  It 


measured  average  and  maximum  tin*  in  iha  system  aa  wall  aa  raaourca 
utilization.  The  system  required  approximately  10  CRAF  aircraft  and 
between  14  and  17  C-9  aircraft.  Tha  average  time  in  tha  system  waa 
undar  thraa  and  ona-half  houra. 

Bacauaa  of  tha  diffaranca  batwaan  tha  praaant  CONUS-based  Air  Forca 
C-9  invantory  and  tha  number  raquirad  by  tha  simulation,  thay 
recommended  altarnata  systems  ba  davalopad.  Ona  of  thaaa  includad  tha 
axcluaiva  uaa  of  CRAF  and  C-9  aircraft  for  avacuation,  tha  basia  of  thia 
raaaarch  atudy. 

Two  othar  raaaarch  afforta  employed  simulation  tachniquaa  to  modal 
a  complex  avacuation  problem.  Honcura  and  Uhita  (24)  and  Gullatt  and 
Stivar  (14)  invaatigatad  tha  Non-combatant  Evacuation  Oparation  (NEO) 
for  tha  Fadaral  Rapublic  of  Germany  in  casa  of  a  major  conflict. 

Tha  NEO  plan  antaila  tha  avacuation  of  approximataly  750,000 
military  dependents,  atata  dapartmant  employees,  and  othara  from  Uaat 
Garmany.  Tha  atudiaa  convartad  all  CRAF  and  C-5  aircraft  to  C-141 
aquivalant  aircraft.  Than,  for  varioua  capacitiaa  and  aircraft 
availabilitiaa,  tha  modala  providad  avacuation  timaa  that  wara  compared 
to  area  "overrun11  times  under  different  scenarios.  Both  research 
efforts  used  Q~Gert  simulation  techniques  to  model  the  system. 

>  As  in  the  previous  studies,  discrete  event  simulation  techniques 

will  be  used  in  this  study  to  modal  tha  network.  Tha  state  of  the 
system  changes  only  with  the  occurrence  of  specific  events.  Some  of 
these  events  include  the  release  of  patients  from  hospitals  in  Europe, 
their  arrival  in  the  U.S.  and  their  discharge  from  the  recovery 
hospital.  Each  significant  event  and  intervening  activity  duration  will 
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be  sinulated  in  th*  Model  using  Sinulntion  Language  for  Alternative 
Modeling  (SLAM)  by  Pritsker  (3012). 

Sctoacifi 

Th*  strategic  evacuation  and  CONUS  distribution  of  casualties 
resulting  froM  a  conventional  war  in  Europe  will  b*  sinulated  for  60 
days.  This  is  th*  san*  tin*  p*riod  used  by  Alfano  and  O'Neill  and  by 
th*  Arn*d  Services  H*dical  Readiness  Office  for  nodical  planning  and 
fore*  requirenents  (8).  Th*  scenario  assun*s  that  by  th*  sixtieth  day, 
s*a  linos  of  r*supply  will  hav*  been  established  thus  reducing  the 
airlift  requirenent.  These  airlift  resources  could  then  be  used  for 
evacuation. 

During  the  first  60-day  period,  this  study  assunes  that  C-141 
aircraft  will  not  be  used  for  casualty  evacuation.  CRAF  aircraft, 
prinarily  6-767,  will  evacuate  the  casualties  and  augnent  the  C-9 
aircraft  in  the  distribution  phase.  If  the  C-9  aircraft  are  deployed, 
other  nodified  CRAF  or  nore  B-767  aircraft  would  be  required. 

Meiwsck  Sicuciuce 

Qolfild*  The  original  Casualty  Evacuation  Model  (CEH)  had  nine 
onload  locations  throughout  Europe.  In  this  study,  these  locations  have 
been  reduced  to  one  centrally  located  point.  The  prinary  Measure  of 
effectiveness  (average  tiMe  in  the  systen )  is  not  affected  by  this 
amplification  because  the  patient's  tiMe  does  not  start  until  a  full 
load  of  patients  is  released  for  evacuation.  The  amplification 
eliiiinates  the  need  to  specify  the  nine  onload  points  and  the  proportion 
of  total  casualties  departing  froM  each  onload  point. 

18 

■"iMn k  nvn.k *  >  •  wtvp  .  «.k  *v» v  ■  *■»  V*  •  /> »• 


UliMMT ENJt-V  -V 


CQEfiUt*  Iiu*.  The  corresponding  flying  tiMes  froH  nine  locations 
has  also  been  eliMinated.  Sanples  fron  ona  of  four  uniforn 
distributions  with  paranatars  basad  on  avaraga  Mininun  and  MaxinuM 
flying  tiMes  to  aach  of  four  U.S.  staging  basas  will  raprasant  tha 
flying  tiMes  for  aach  transatlantic  Mission. 

EnCfiUlt  Slags.  Raplacing  tha  original  74  intarnadiata  locations  of 
tha  original  Modal  (20)  ara  tha  four  staging  basas  -  Boston,  Naw  York, 
Philadelphia  and  Washington.  In  tha  original  CEH,  avary  likaly  anrouta 
stop  was  includad  bacausa  of  tha  Many  cargo  onload  locations  and 
nuMarous  routings.  Bacausa  tha  CRAF  Mission  is  exclusively  a  nodical 
Mission  in  this  study,  thasa  anrouta  locations  ara  not  required. 

Tha  spacific  staging  basas  wara  salactad  primarily  dua  to  tha  larga 
nuMbar  of  hospitals  naarby.  Thay  ara  also  within  tha  B-767  fual  ranga. 

Each  location  prasantly  supports  intarnational  flights  and  is  a  Major 
oparations  cantar  for  Many  coMMarcial  carriers.  Activation  of  tha  CRAF 
plan  and  tha  nodical  tasking  would  impose  Mininal  ground  support  changas 
to  tha  axisting  structure. 

C"2  QptCltifiDS*  Interface  between  tha  B-767  and  C-9  aircraft  will 
occur  only  in  Washington.  Tha  C-9  distribution  network  includes  all 
hospitals  located  within  about  two  hours  flying  tine  ( approxinately  900 
Milas )  fron  Washington  except  those  hosoitals  served  by  tha  cthe*'  three 
staging  bases.  Tha  seni-circle  in  Figure  3  roughly  outlines  the  C-9 
distribution  area. 

B~ZfiZ  QptESlifiDl.  Boston,  Naw  York  and  Philadalphia  will  serve 
hospitals  located  within  approxinately  2  hours  driving  tine  (about  100 
Milas).  Figure  3  identifies  the  areas  in  the  western  U.S.  served  by 
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each  staging  bats.  The  dashtd  line  indicatas  tha  extent  of  8-767 
oparations  required  for  this  nodal. 

Data  Eilt-  An  array  of  citias  served  by  aach  staging  basa  is 
ordarad  according  to  its  distanca  fron  tha  stags  location.  Tha  array 
data  includas  tha  flying  tins  and  nunber  of  availabla  hospital  bads  for 
aach  city.  Tha  data  is  part  of  tha  initialization  subroutina  of  tha 
fortran  conputer  coda  shown  in  Appandix  A. 


Cities  can  ba  addad  or  dalatad  fron  tha  network  with  only  ninor 
changas  to  tha  coda,  but  tha  ordarad  saquanca  basad  on  distanca  nust  ba 
naintainad.  This  will  ansura  that  tha  closast  citias  with  sufficiant 
anpty  bads  ara  salactad  first  for  aach  staging  basa  by  tha  fortran 
subroutina  "3ERCH". 


CB6E  HesiiaaliQD-  Tha  proportion  of  B-767  Missions  destined  to 
aach  staging  basa  was  datarninad  fron  tha  total  nunber  of  bads  sarvad  by 
that  basa  and  tha  nunber  of  bads  within  two  hours  driving  tine.  Tha 
requirenent  was  that  aach  area  ba  filled  at  approxinataly  tha  sane  rata. 


Half  of  all  CRAF  nissions  ara  scheduled  to  Washington;  tha  other 
half  ara  divided  anong  tha  other  three  citias.  This  avoids  saturation 
of  ground  facilities  and  provides  sona  separation  between  arriving 
flights.  Every  other  B-767  fron  Europe  is  directed  to  Washington  by  tha 
"SELECT"  node  in  tha  SLAM  network  code. 


Whan  an  aircraft  arrives  at  Boston,  New  York  or  Philadelphia  and 
tha  nunber  of  anpty  bads  in  tha  innadiate  area  is  equal  to  or  greater 
tha  capacity  of  the  B-767,  all  patients  are  offloaded  and  the  B-767  is 
returned  to  Europe.  Otherwise,  the  Mission  continues  to  tha  first  city 
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served  by  th*  tinging  bat*  with  enough  enpty  bed*  for  all  patient* 
aboard.  Th*  patient*  are  offloaded  and  th*  aircraft  it  returned  to  the 
staging  bate  and  then  to  Europe. 

When  an  aircraft  arrive*  at  Washington,  all  casualties  are 
offloaded  and  th*  CRAF  is  returned  to  Europe.  For  each  40  patients 
offloaded  (C-9  capacity),  two  conditions  are  checked.  For  th*  first 
condition,  a  subroutine  searches  an  array  for  th*  nearest  city  with  at 
least  40  available  beds.  Th*  second  condition  requires  at  least  one  C-9 
be  available.  If  both  conditions  are  net,  th*  patients  are  flown  to  th* 
selected  city. 

ElilftOi  Qk*c£lfiW*  For  all  flights  to  the  U.S.,  an  overflow 
condition  can  develop  when  on*  of  th*  following  situations  occurs! 

1.  If  no  beds  are  available  in  Boston,  New  York  or  Philadelphia 
area  hospitals  or  their  respective  service  areas,  or 

2.  If  no  C-9  aircraft  are  available  at  Washington  and  no  beds  are 
available  in  Washington  area  hospitals,  or 

3.  If  no  beds  are  available  in  the  C-9  servicing  area  or  in  th* 
Washington  area  hospitals,  or 

4.  If  less  than  30  patients  renain  at  Washington  after  full  C-9 
loads  have  departed  and  no  beds  are  available  in  Washington 
area  hospitals. 

In  these  cases,  all  patients  are  offloaded,  adnitted  to  area 
hospitals  and  assigned  to  beds  that  were  not  previously  contracted  for 
use  through  th*  National  Disaster  Hedical  Systen  (NMDS).  Recall  that 
25%  of  the  beds  in  each  participating  hospital  are  contracted.  If  the 
overflow  patients  require  an  additional  10%  of  a  hospital's  capability 
( total  35% ),  then  th*  selected  conbination  of  input  paraneters  is 
unacceptable  and  will  be  rejected  as  a  viable  option. 
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Th*  nodal  Maintain*  statistic*  on  th#  nunber  of  patients  causing  an 
ovarflow  conditions  and  where  it  occurs.  In  tasting  tha  nodal  at  high 
casualty  ratas,  ovarflow  occurrad  at  Washington  but  tha  nunber  of 
casual tias  causing  this  condition  was  lass  than  tha  Units  shown  in 
Tabla  I  basad  on  tha  35Z  rastriction  abova. 


Tabla  I 

Hospital  Ovarflow  Data 


Stag*  Location 

Contracted  Bads 

Ovarflow 

Washington 

4360 

1740 

New  York 

13420 

5360 

Boston 

4410 

1760 

Philadelphia 

4390 

1750 

Ealital  CstftQocita 

Tha  Arnad  Services  Hadical  Raadinass  Offica  has  classifiad  all 
conbat  casualtias  into  309  catagorias  (8>.  Tha  Hilitary  Airlift  Connand 
has  aggragatad  thasa  catagorias  into  tha  classifications  shown  in  Tabla j 
II  (U23).  Tha  astinatad  parcantaga  of  total  casualtias  in  aach 

i 

catagory  and  tha  expected  nunbar  of  day*  in  racovary  is  included.  In 
casa  tha  racovary  periods  have  bean  underestinated,  a  252  increase  in 
tha  racovary  period  will  be  tasted  during  tha  nodal  exparinantation 
phase.  Tha  results  will  be  included  in  chapter  IV. 

Thasa  patient  catagorias  have  bean  reduced  based  on  their  length  of 
stay  in  th*  racovary  hospital  for  this  study.  It  assunes  that  every 
hospital  can  adnit  any  catagory  patient?  however,  son*  casualty 
screening  process  could  be  inplanented  for  certain  cases.  Th*  revised 
categories  ar#  shown  in  Table  III. 


Assuning  every  hospital  can  adnit  any  patiant  elininates  th*  need 


to  Match  a  patiant  with  a  particualar  typa  hospital  bad.  Tha  nodal 
determines  discharga  datas  basad  on  patiant  catagory  and  tha  arrival 
data. 


Tabla  II 

Patiant  Charactaristics 


Catagory  X  of  Total 

Recovery  Period 

1.  Hadical 

20 

1A 

2.  Psychiatric 

7 

29 

3.  General  Surgery 

31 

24 

4.  Orthopedic  Surgery 

19 

50 

3.  Neurosurgery 

6 

36 

4.  Oral/Haxillo  Facial 

7 

40 

7.  Urology 

1 

1.2 

Q.  Opthanology 

3 

27 

7.  Burns 

2 

33 

10.  Thoracic  Surgery 

4 

54 

11.  Spinal  Cord 

(<  1 

38 

Tabla  III 


Revised  Patiant  Charactaristics 


New 

Old 

Recovery 

Catagory 

Catagory 

Percent 

Period 

I 

1.  7 

21 

16 

II 

2,  3f  8 

41 

23 

III 

5,  4,  9,  11 

15 

38 

IV 

4,  10 

23 

51 

daautiBii&DS  iad  libit a ticca 

Assunptions  and  linitations  ara  nacassary  to  sinplify  tha  nodal  and 
raduca  tha  nunber  of  computations.  Sone  apply  to  tha  structure  of  the 
nodal  while  others  refer  to  tha  data.  While  sone  have  already  been 
Mentioned  and  others  have  not  yet  bean  discussed,  all  assunptions  and 
linitations  will  be  listed  here  for  easier  reference. 


Tactical  Medical  avacuation  operations  are  presumed  functioning 
properly;  therefore,  CRAF  aircraft  will  not  be  subject  to  the  hazards 
of  flying  near  hostile  areas. 

HlinltOICCt*  Mechanical  delays  and  diversions  are  not  Modelled  in 
the  simulation;  however,  additional  ground  tine  is  scheduled  for  each 
eastbound  CRAF  at  the  stage  location.  An  adjustiient  in  the  mmber  of 
aircraft  required  by  the  Model  will  be  proposed  in  the  next  chapter  to 
account  for  the  absence  of  Maintenance  factors.  An  increase  in  the 
nunber  of  B -767  aircraft  in  tho  systeM  will  offset  the  increase  in  the 
patient's  average  tine  in  the  systeM  caused  by  Maintenance  delays  and 
diversions. 

For  the  C-?  aircraft,  no  More  than  nine  are  used  in  the  Model  even 
though  the  Air  Force  has  11  located  in  the  li.S.  The  extra  aircraft  can 
be  used  for  Maintenance  spares  or  for  special  Medical  Missions  not 
included  in  this  study. 

licatbcc.  Belays  and  diversions  caused  by  weather  have  not  been 
incorporated  into  the  Model.  These  deviations  would  have  to  be  frequent 
and  extensive  to  Make  a  statistically  significant  difference  in  the 
patient's  average  tine  in  the  systeM.  A  Modification  siMilar  to  the 
Maintence  adjustment  will  be  proposed  to  allow  for  this  factor. 

Aicccafii  Utilization.  The  Model  requires  that  all  B-767  aircraft 
be  fully  loaded  before  take-off.  This  MaxiMizes  aircraft  capacity 
utilization  without  affecting  the  primary  Measure  of  systeM  performance 
or  causing  any  additional  delay  of  patients  already  waiting  for 
transportation.  Also,  continuation  flights  for  B-767  aircraft  arriving 
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at  Boston ,  Now  York  or  Philadelphia,  will  bo  H«do  only  to  destination 
cities  that  can  adnit  a  full  load  of  patients  for  Maxinun  aircraft 
efficiency. 

The  Model  does  allow  the  C-9  to  depart  with  75X  of  its  nornal 
capacity  in  order  to  reduce  the  lumber  of  overflow  patients  at 
Washington  area  hospitals.  Thus,  if  110  casualties  arrive  on  one  B -767, 
30  patients  would  reMain  after  two  full  C-9  loads  have  departed.  These 
patients  would  fill  Washington  area  hospitals  within  20  days  if  eight 

I 

aircraft  arrived  each  day.  If  130  casualties  arrive  on  each  B-767,  only 
10  patients  would  renain  after  three  full  C-9  loads  have  departed.  For 
the  sane  arrival  rate  as  above,  Washington  area  hospitals  would  not  be 
filled  within  60  days. 

fiCfiUDd  Supped.  Ground  transportation  required  to  Move  patients 
between  hospitals  and  aircraft  or  between  aircraft  is  assuwed  available. 
This  includes  vehicles  for  the  trips  to  hospitals  within  approxinately 
100  Miles  of  any  offload  location. 

CBBE  Capacity*  The  passenger  capacity  of  a  B-767  is  about  220 
(34)|  however,  a  litter  patient  occupies  about  two  and  one-half  tines 
More  space  than  an  anbulatory  patient.  The  expected  nunber  of  litter 
patients  varies  fron  602  to  902  of  the  total  nunber  of  casualties.  This 
reduces  the  patient  capacity  to  between  100  and  140  for  the  CRAF.  These 
factor  levels  will  be  discussed  in  the  next  chapter. 

Caauallitl.  Estinates  of  the  nunber  of  casualties  per  day  are 
classified)  however,  Alfano  and  O'Neill's  study  and  other  unclassified 
estinates  have  used  1000  or  More  per  day  (8).  To  Model  a  variable 
casualty  rate,  a  triangular  distribution  with  a  low  paraneter  value  of 
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1000  and  a  high  paraneter  valua  of  2500  will  ba  u*ad  to  provide  a 
raaaonably  wida  range  of  casualty  rata*.  Tha  triangular  distribution 
was  salactad  because  of  its  sinplicity.  Also,  it  does  not  have  tha 
problen  associated  with  extrene  values  found  in  tha  tails  of  other 
distributions.  Tha  affect  of  various  nodal  values  will  ba  exanined  in 
tha  next  chapter. 

UfilBiiSll*  All  hospitals  located  within  two  hours  driving  tine 
(about  100  nilas)  fron  each  offload  location  have  bean  aggregated  in  tha 
array  of  cities  for  each  staging  base. 

CfitJPUttC  line.  To  save  conputar  processing  tine,  tha  nodal  will 
generate  one  entity  for  every  tan  casualties.  Hospital  bad  availability 
and  aircraft  capacities  have  bean  appropriately  reduced.  Thus,  tha 
capacity  of  a  C-?  will  ba  4  instead  of  40  and  a  B -767  aircraft  capacity 
will  range  fron  10  to  14. 

ytcificatiQD 

Hodel  verification  is  an  inherent  part  of  nodal  developnent  fron 
its  conceptualization  through  sunnary  data  analysis.  Verification 
ensures  that  the  conceptual  nodal  is  accurately  reflected  in  the  SLAH 
and  fortran  conputar  code  (4*379).  One  other  result  of  verification  is 
increased  confidence  in  the  nodal. 

The  following  verification  techniques  were  carried  out* 

1.  Logic  verification  -  the  conputer  code  was  critiqued  by  ny 
thesis  advisor  for  logic  and  progranning  nethods. 

2.  Flow  diagran  -  entities  created  in  the  network  were  traced 
through  the  systen.  Every  possible  path  terninated  at  a  data 
collection  point  so  no  entity  can  be  trapped  in  the  network. 

3.  Docunentation  -  the  nain  progran  and  each  nodule  of  the 
conputer  code  include  a  general  description  of  its  purpose. 

Each  new  variable  is  defined. 
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4.  Sensitivity  testing  -  each  paraneter,  in  turn,  was  varied  while 
all  other i  were  fixed  to  verify  that  the  behavior  of  the  nodel 
was  at  expected.  The  concern  here  was  to  verify  that  the 
response  variable  increased  or  decreased  when  it  should  without 
regard  for  the  anount  of  increase  or  decrease. 

5.  Stress  testing  -  each  paraneter  was  independently  tested  at 
high  and  low  values  to  cause  the  nodel  to  "blow-up"  or 
stabilize  as  expected. 


UalidatifiQ 

No  nodel  is  ever  totally  representative  of  the  systen  (4*384),  but 
the  ultinate  test  of  its  validity  is  its  ability  to  predict  the  future 
behavior  of  the  real-world  systen.  For  the  scenario  sinulated  in  this 
study,  it  will  hopefully  never  happen)  however,  other  validating 
techniques  are  available  to  increase  the  Modeler's  confidence  that  the 
nodel  is  accurate.  The  following  techniques  were  applied* 


1.  Face  validity  -  the  nodel  appears  reasonable.  It  behaves  as 
expected  when  input  variables  are  changed.  A  high  degree  of 
realisn  is  built  into  the  nodel  through  reasonable  assunptions 
regarding  the  network  structure  and  reliable  data. 

2.  Structural  assunptions  -  the  Model's  sinplifications  and 
abstractions  of  the  evacuation  process  appear  reasonable  based 
on  ny  prior  experience  in  strategic  and  tactical  airlift 
operations. 

3.  Sensitivity  analysis  -  the  nodel  provides  reasonable  data  when 
input  paraneters  are  varied.  Unlike  sensitivity  testing  for 
verification,  sensitivity  analysis  is  concerned  with  the  anount 
or  degree  of  change  in  the  response  variable. 


Suatttty 

Sinulation  techniques  are  ideally  suited  to  evaluate  the 
perfornance  of  an  evacuation  and  distribution  network.  Fortran  conputer 
code  and  Sinulation  Language  for  Alternative  Modeling  (SLAM)  will  be 
conbined  in  a  network-discrete  event  sinulation  to  nodel  the  evacuation 
and  distribution  of  casualties  resulting  fron  a  60-day  conventional 
conflict  in  Europe. 
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The  structure  of  the  network  includes  one  onload  location,  four 
staging  bases  and  a  total  of  49  offload  points  with  a  coMbined  total  of 
83180  available  hospital  beds  in  87  cities.  Each  staging  base  serves 
cities  listed  in  an  array  based  on  shortest  distance  first  and  includes 
the  flying  tine  to  each  city  and  the  nunber  of  beds  available  in  its 
hospitals. 

The  original  309  categories  of  patients  has  been  reduced  to  four 
categories  based  on  their  recovery  tines.  Also,  B-767  aircraft  will  be 
fully  loaded  for  every  flight)  no  aircraft  will  proceed  to  destinations 
with  insufficient  beds  for  the  entire  load. 

All  najor  assuMptions  and  linitations  to  reduce  the  conplexity  of 
the  Model  and  coHputations  required  have  been  included.  These  refer  to 
both  structural  and  data  assuMptions,  sone  of  which  will  be  further 
explained  in  the  experinental  design  of  the  Model  in  the  next  chapter. 


Ill.  EXPERIMENTAL  DESIGN 

iD&EfidUCiiQD 

The  Casualty  Evacuation  Nodal  ( CEH )  davalopad  in  chapter  II 
similates  a  systeri  conposed  of  factors  that  affact  a  patient's  avarage 
tin*  in  tha  evacuation  network.  This  chaptar  will  davalop  an 
axparinantal  dati$.'  us  des.ribed  by  Montgonary  (25)  using  analysis  of 
varianca  ( ANOVA )  tachriquas  <  16  )  to  quantify  tha  statistical 
significanca  of  tha  factors  and  thair  intaractions. 

Tha  purpose  of  tha  ANOVA  is  to  coMpare  tha  diffarant  system 
spacifiad  by  tha  axparinantal  design  and  determine  any  factors  or 
intaractions  that  can  ba  eliminated  fron  tha  nodal  or  fixad  at  tone 
arbitrary,  Lut  reasonable,  level*.  A  factor  nay  have  no  statistical  or 
practical  significance  ir  tha  *  odal  if  a  change  in  tha  level  of  tha 
factor  causes  no  apore-r>.  change  in  tha  response  variable  (average 
tine  in  th«  system),  when  all  other  factors  are  held  constant. 

After  the  significant  factors  and  interactions  have  been 
datarninad,  Multiple  regression  techniques  will  ba  used  to  davalop  a 
response  surface  equation.  This  will  ba  addressed  separately  in  Chapter 
IV.  The  equation  can  be  used  to  predict  the  value  of  tha  patient's  tine 
in  tha  systeM  given  any  conbination  of  factor  settings  within  a 
specified  range. 

Tha  experiMental  design  developnent  begins  with  data  collection  and 
includes  a  factor  screening  process.  Specific  levels  of  each  factor 
will  be  established  and  coded  depending  on  the  experiMental  design 
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selected.  Other  areas,  including  variance  reduction  techniques  will  be 

addressed. 

Dili  CfilltciifiO 

Chapter  II  included  data  for  the  casualty  rate  distribution  and 
hospitals  that  will  be  used  in  the  simulation  Model.  This  data,  as  well 
as  other  Factor  data,  will  be  included  in  this  section. 

Hfifpiial  LfiCltifiOS*  The  HAC  Operations  Research  Division  provided 
a  list  of  ?4  Metropolitan  areas  and  the  total  nuMber  oF  beds  in  each 
area  as  oF  Mid-l?0S  (?)  shown  in  Appendix  B.  The  locations  were 
aggregated  into  approxinately  40  areas  that  will  be  served  by  the 
staging  bases  outlined  in  chapter  II. 

CfitltJiiitd  Btdl'  The  list  above  included  over  400,000  available 
beds;  however,  only  about  252  oF  these  beds  can  be  coMMitted  For  use  by 
the  Military  (2?).  This  provides  100,000  beds  but  only  85,000  beds 
located  in  the  eastern  halF  oF  the  country  are  necessary  For  this  study. 

Caiualty  Bits-  As  stated  in  chapter  II,  a  triangular  distribution 
with  low  and  high  paraMeter  values  of  1000  and  2500  respectively,  will 
be  used  to  provide  a  wide  range  oF  casualty  rates.  By  changing  the 
Modal  value  of  the  distribution,  the  expected  value  can  be  increased  or 
decreased  as  necessary.  Specific  values  of  the  node  will  be  discussed 
in  a  later  section. 

C86E  Capacity.  Also  fron  chapter  II,  the  passenger  capacity  for  a 
B -767  aircraft  is  about  220;  but  the  space  required  for  a  litter  patient 
is  about  two  and  one-half  tines  that  required  for  an  anbulatory  patient . 
Estinates  of  the  nuMber  of  litter  patients  vary  fron  60X  to  90X  of  the 
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total  casualties.  Thai*  percentages  correspond  to  aircraft  patiant 
capacitiaa  batwaan  100  and  140  based  on  tha  220  passenger  capacity  of 


tha  aircraft.  Thasa  values  rapraaant  high  and  low  linitsj  however, 
othar  capacitiaa  will  ba  avaluatad. 

MuobtE  of  B~Z6Z*  Thia  factor  ia  vary  dependent  on  tha  caaualty 
rata  and  tha  capacity  of  tha  aircraft.  U&ing  raaaonabla  low  and  high 
valuaa  for  tha  caaualty  rata  and  an  avaraga  120  patiant  capacity,  tha 
expected  nunber  of  aircraft  raquirad  ia  batwaan  14  and  16.  Baaad  on 
initial  pilot  runa  of  tha  ainulation  nodal,  highar  valuaa  wara  raquirad. 
This  will  alao  ba  addraaaad  in  a  latar  aaction. 

MuabtC  of  C~2*  Bacauaa  tha  Air  Forca  only  haa  11  C-9  aircraft 
atationad  in  tha  CONUS  (10),  thia  will  ba  tha  naxinun  nunbar  available. 
Howavar,  tha  nodal  will  ba  axarciaad  at  lowar  valuaa  to  datarnina  the 
effect  thia  has  on  tha  rasponaa  variable. 

IcaDiatlaotiC  Elyiog  line.  Transatlantic  crossing  tines  vary  with 
the  route  selected,  altitude  and  tine  of  year.  Banks  and  Carson 
reconnend  using  a  Uniforn  distribution  whan  vary  little  is  known  about 
tha  actual  distribution  or  when  tha  events  are  truly  randon  ( 4 ).  Tha 
paraneters  for  tha  distributions  will  ba  tha  average  nininun  and  avaraga 
naxinun  flying  tines  based  on  tha  aircraft's  cruise  speed  and  average 
wind  factors. 

The  distribution  data  for  the  transatlantic  flights  to  the  four 
arrival  or  staging  bases  is  shown  in  Table  IV.  The  return  flight  to 
Europe  has  a  Uniforn  distribution  also  shown  in  Table  IV. 
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CQtJUS  Eli'iog  list.  The  flying  tines  for  all  CONUS  flights  are 
based  on  the  distance  between  airports,  aircraft  speed  and  wind  factor} 
but  the  effect  of  short  distances,  where  the  aircraft  never  attains  its 
design  cruise  speed,  is  also  incorporated.  This  data  is  stored  in  the 
array  of  cities  served  by  each  staging  base  shown  in  Appendix  A. 

Table  IV 

B-767  Flying  Tine  Distributions 


Destination 

Flying  Time 

Boston 

Unfrm( 7.7,  B.7> 

New  York 

Unfrm( 8.0,  9.0) 

Philadelphia 

Unfrn<  8.2,  9.2) 

Washington 

Unfrn<  8.5,  9.5) 

Europe 

Unfrm(7.0,  8.5) 

Gcfiuod  lias.  Nornal  ground  tine  for  onload,  offload  and  refueling 
during  contingency  operations  is  two  and  one-half  hours  for  C-141 
aeronedical  nissions.  No  guidelines  are  available  for  nissions  with  as 
nany  as,  say  120  patients,  but  reasonable  estinates  can  be  nade  using 
the  nornal  C-141  tine  assuning  a  full  load  of  67  patients.  The  average 
tine  per  patient  is  about  two  ninutes  and  13  seconds.  For  the  B-767 
with  120  patients  aboard  (732  litter),  this  rate  would  require  3.375 
hours.  C-?  ground  tine  can  be  deternined  similarly. 

Assuning  all  ground  operations  can  be  conducted  concurrently  with 
onloading  or  offloading,  tho  data  in  Table  V  will  be  used  in  the 
simulation  model.  Additional  tine  is  provided  for  eastbound  CRAF  at 
staging  locations  for  minor  maintenance.  Purpose  codes  are* 

on  ■  onload  cc  *  crew  change 

off  ■  offload  ret'  *  refueling 

nx  *  maintenance 
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Table  V 


B -767  and  C-?  Ground  Tinas  (Hours) 


«  * 
i 


Aircraft 

Departure 

Tine 

Purpose 

CRAF 

Europe 

3.5 

On,  ref,  cc. 

II 

West  Stage 

1.5 

Ref,  cc. 

It 

East  Stage 

4.0 

Off,  ref,  cc,  nx. 

It 

Other  points 

3.0 

Off,  ref. 

C-? 

Washington 

2.5 

On,  ref,  cc. 

It 

Other  points 

2.0 

Off,  ref. 

Initial  Scctsoisfi 

All  factors  discussed  above  are  inportant  aspects  of  the  real-world 
systen  sinulated  in  the  nodelf  however,  one  or  nore  of  then  nay  not  be 
significant  in  terns  of  changes  in  the  response  variable.  Other  factors 
directly  affect  the  patient's  tine  in  the  systen,  but  are  generally 
beyond  our  control  or  should  renain  at  fixed  levels.  The  specific 
factors  in  these  case  are  the  flying  tine  and  ground  tine. 

For  the  experinental  design,  the  flyinQ  tine  and  ground  tine 
factors  will  not  be  evaluated  at  high  or  low  values  to  exanine  their 
effects  on  the  response  variable.  The  CONUS  flying  tine  will  renain 
constant  while  the  transatlantic  tines  will  be  selected  fron  the  listed 
distributions. 


The  ground  tines  will  also  renain  fixed  in  order  to  allow  an 
orderly  conpletion  of  all  ground  activities  in  the  proper  sequence. 
While  shorter  ground  tines  reduce  a  patient's  tine  in  the  systen,  the 
intense  activities  required  by  Haintenance,  aircrews,  nedcrews,  aerial 
port  teans  and  conHand  post  controllers  to  attain  shorter  ground  tines 
could  not  be  sustained  for  a  very  long  tine. 


; 


-  By  including  flying  tin*  and  ground  tin*  a*  fixed  factors  within 

th*  sinulation  nodal,  four  factors  renain  in  th*  experinental  design. 
High  and  low  values  of  *ach  factor  will  b*  specified  next. 

Eactfic  Ltyala 

Casualty  Bit*.  A  wide  rang*  of  daily  casualty  rat*s  is  guaranteed 
by  th*  choic*  of  th*  high  and  low  paran*t*r  values  of  th*  distribution. 
By  setting  th*  nodal  valu*  at  relatively  high  or  low  values,  the 
expected  values  will  be  high  or  low  accordingly.  For  th*  experinental 
design,  a  nodal  value  of  1300  will  result  in  an  expected  value  of  1600 
casualties  per  day.  A  nodal  valu*  of  2200  will  cause  th*  expected  valu* 
to  reach  1700  per  day.  In  both  cases,  a  sanpl*  fron  th*  distribution 

j 

I  can  still  have  very  high  or  very  low  values,  depending  on  the  randon 
nunber  selected. 

CBAE  Capacity.  As  stated  earlier,  th*  proportion  of  litter 
patients  anong  all  casualties  ranges  fron  60%  to  90%.  The  HAC  Surgeon 
l  General's  Office  has  used  a  patient  capacity  of  about  120  for  son* 

nodical  planning  studies  (27).  This  suggests  that  a  75%  litter  rate  and 
a  B-767  passenger  capacity  of  220  seen  reasonable.  So,  to  deternine  if 
th*  capacity  of  th*  aircraft  affects  the  response  variable  or  if 
interactions  are  significant,  th*  experinental  design  will  be  conducted 
with  aircraft  capacities  of  110  and  130. 

ttUBbtC  s£  B~Z£Z*  Based  on  pilot  runs  of  the  sinulation  nodei,  a 
reasonable  range  of  the  nunber  of  aircraft  to  be  considered  in  the 
design  was  fron  16  to  18.  The  original  range  of  14  to  16  resulted  in 
high  response  variable  values  caused  by  too  few  aircraft  and  too  nany 
patients  waiting  for  transportation. 
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NuubiC  o£  G*2*  The  design  will  exanine  the  effects  of  using  only 
six  or  nine  C-?  aircraft.  This  will  also  provide  an  indication  of  how 
Many  aircraft  would  be  required  to  replace  the  C-?'s  should  they  be 
deployed  to  Europe. 

EiDfll  Otiigs 

A  full  factorial  design  will  be  developed  to  quantify  the  Main 
effects  of  each  factor  and  the  inter  ictions  between  factors.  Based  on 
the  pilot  runs  of  the  Model  and  the  initial  screening  process,  the  low 
and  high  levels  of  each  factor  have  been  deternined  and  are  sumiarized 
in  Table  VI. 


Table  VI 

Experinental  Design  Factor  Settings 


NQDzCQnstaDi  Variance 

The  Model  sixulates  a  potential  queueing  situation  under  certain 
treatMent  settings.  In  such  cases,  the  response  variable  will  be 
significantly  autocorrelated  (16:438).  The  reason  is  clear:  the  tine  in 
the  systex  for  an  arriving  patient  depends  on  the  tine  in  the  systen  for 
a  patients  arriving  ahead  of  hin.  The  Model  queue  discipline  is  a 
first-in  first-out  process. 

The  detriMental  effects  of  autocorrelated  data  can  be  reduced  by 
appropriate  data  transforxations  as  suggested  by  Draper  and  SMith  ( 11 ) 
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and  Montgomery  (25).  A  well-known  transformation  and  the  one  used  in 
this  model  is  the  natural  logarithm  transformation  of  the  response 
variable.  The  resulting  model  is  an  exponential  model  in  the  original 
response  variable.  The  logarithm  of  the  random  errors  have  a  Normal 
distribution  as  a  result  of  the  transformation. 

yaciaoct  BeduclieD 

One  simple  variance  reduction  technique  is  blocking  on  common 
random  number  streams  (41458).  Each  distribution  in  the  model  will  use 
a  unique  random  number  generator  provided  by  SLAM.  As  a  result, 
observed  differences  between  means  will  be  due  to  actual  differences 
rather  than  due  to  the  variance  associated  with  the  distributions. 

Variance  reduction  can  also  be  achieved  by  use  of  antithetic 
variates  by  inducing  negative  correlation  in  the  response  data  (30i506). 
Unfortunately,  two  simulation  runs  are  required  to  produce  a  single 
observation.  Also,  Law  and  Kelton  do  not  recommend  mixing  the  common 
random  number  stream  technique  with  the  antithetic  variate  technique, 
since  cross  covariances  might  actually  increase  the  variance  (18).  As  a 
result,  antithetic  variates  will  not  be  used  in  this  simulation  study. 

Sauple  Size 

Montgomery  provides  a  procedure  fof  estimating  main  effects  and 
interactions  when  only  a  single  replicate  of  each  treatment  can  be  made 
(25J273).  The  method  assumes  that  certain  high-order  interactions  are 
negligible  and  may  be  combined  to  estimate  the  experimental  error. 
Because  there  are  only  16  treatment  combinations  in  this  design,  two 
replicates  were  made  for  each  treatment  in  order  to  obtain  a  more 
accurate  estimate  of  the  experimental  error. 


Analysis  q£  UaciaDce 

Because  only  two  levels  of  aach  factor  ara  evaluated,  tha  response 
variabla  it  atiuMad  to  ba  linaar  ovar  tha  range  of  tha  factor  lavalt 
aalactad.  However,  tha  raaponaa  variabla  for  tha  experimental  design  it 
tha  logarithm  of  tha  simulation  raaponaa  value;  tharafora  tha  logarithm 
of  tha  raaponaa  variabla  ia  assumed  to  ba  linaar.  It  may  ba  inaccurata 
in  dascribing  what  ia  assantially  a  queueing  procass,  but  it  will  ba 
aufficiantly  cloaa  in  ordar  to  datarmina  significant  factors  and 
intaractiona. 

Conlcsst  Constants.  Tha  tablas  of  experimental  dasign  algabraic 
sign  for  calculating  tha  factorial  affects  for  tha  3-factor  and  4-factor 
modala  can  ba  found  in  Montgomery  (  25 J268, 277  ).  A  statistical  package 
( BMDP )  was  used  with  these  contrast  constants  and  tha  transformed 
response  variables  to  datarmina  the  factorial  effects. 

Tha  four -factor  analysis  of  variance  table  is  shown  in  Appendix  C. 
Significant  factors  and  interactions  are  clearly  evident  from  the  F- 
statistic.  The  number  of  C-9  aircraft  and  all  interactions  involving 
that  factor  are  statistically  insignificant  in  terms  of  changes  in  the 
response  variable.  As  a  result,  the  number  of  C-9  aircraft  in  the 
simulation  model  will  ba  arbitrarily  fixed  at  tha  higher  value. 

An  analysis  of  variance  was  conducted  on  the  three  remaining 
factors  with  an  extra  replicate  included  and  is  also  show  in  Appendix  C. 
All  interactions  remain  significant.  The  results  of  the  experimental 
design  can  be  used  to  develop  a  response  surface  equation  to  be 
discussed  in  the  next  chapter. 
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Exparinantal  design  techniques  enable  a  tiodalar  to  detarnin*  those 
factor*  and  intar action*  that  ar*  statistically  significant  by  conparing 
th#  changas  in  tha  rasponsa  variabla.  Following  tha  data  collaction  and 
nodal  davalopnant  phases,  th*  evacuation  nodal  includad  six  inportant 
factors  that  affact  a  patient's  avaraga  tin*  in  tha  systan. 

An  initial  scraaning  of  thasa  factors  rasultad  in  satting  tha 
flying  tin*  and  ground  tin*  factors  at  spacifiad  distributions  or 
constant  valuas.  Specific  lavals  war*  salactad  for  tha  ranaing  four 
factors  in  tha  final  exparinantal  dasign. 

An  analysis  of  varianca  was  parforxed  on  tha  transforned  rasponsa 
variabla  and  th*  contrast  constants  resulting  in  th*  alinination  of 
another  factor  fron  the  design.  All  three  renaining  factors  -  casualty 
rate,  nunber  of  CRAF  and  CRAF  capacity  -  and  their  interactions  were 
statistically  significant.  Th*  rasponsa  surfaca  aquation  that  can  ba 
davalopad  will  be  discussad  in  Chaptar  IV. 
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XV.  RESPONSE  SURFACES 


XoiCfidUCtifiD 

The  experinental  design  procedures  in  chapter  III  indicated  three 
significant  variables  that  will  be  used  to  deternine  the  patient's  Mean 
tine  in  the  systeM  (TIS)  for  the  evacuation  structure  in  this  study. 
This  chapter  will  use  those  three  variables  -  Mean  casualty  rate,  number 
of  aircraft,  and  aircraft  capacity  -  and  the  statistical  technique 
called  response  surface  Methodology  (RSH)  to  develop  an  equation  that 
will  predict  the  Mean  TIS  value  (28i£2>. 

For  the  casualty  evacuation  Model  ( CEM >,  a  response  surface 
equation  is  valuable  because  it  eliminates  the  need  to  conduct 
additional  siMulation  experinents  in  order  to  deternine  the  patient's 
Mean  TIS  for  any  values  of  the  three  variables  within  established 
bounds.  This  saves  tiMe  and  enables  an  analyst  or  project  officer  to 
exaMine  the  results  of  Many  coMbinations  of  casualty  rates,  nuMber  of 
aircraft  and  aircraft  capacities  to  deterMine  acceptable  options  based 
on  a  subjective  or  liniting  criterion  value  of  the  Mean  TIS. 

The  process  of  generating  a  response  surface  involves  finding  a 
suitable  approxiwation  for  the  true  functional  relationship  between  the 
response  variable  and  the  independent  predictor  variables  (25:445).  The 
functional  forM  of  the  approxinating  polynomial  and  the  estiMation  of 
its  paraMeters  will  oe  discussed  in  this  chapter. 

EuacUaoal  Eqlb 

Because  of  the  transfomation  to  reduce  the  effects  of  nonconstant 
variance  described  in  Chapter  III,  the  forM  of  the  functional 
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relationship  is  exponential  in  th*  original  response  variable;  however, 
the  functional  form  of  the  exponent  nust  be  determined.  After  the 
original  predictor  variables  are  coded  as  in  Chapter  III,  a  linear  model 
will  be  examined. 

Linear:  ticdela 

The  simplest  functional  form  of  the  polynomial  exponent  is  the 
linear  model.  If  the  predictive  accuracy  of  this  form  is  acceptable  to 
the  analyst,  the  exponential  model  with  linear  exponent  can  quickly  and 
easily  provide  useful  estimates  the  the  mean  TIS  for  the  evacuation 
system.  If  greater  precision  is  required,  a  quadratic  model  will  be 
developed. 

Cfidad  Uaniabltl.  Each  independent  variable  can  be  recoded  in  the 
same  manner  used  in  the  experimental  design  techniques  in  chapter  III. 
These  are  the  familiar  contrast  coefficients  coded  from  the  low  and  high 
values  of  each  factor  affecting  the  response  variable  value.  Table  VII 
includes  the  coded  variable  values  and  the  transformed  response  variable 
that  will  be  used  in  developing  the  response  surface  equation. 


Table  VII 


Data 

For  A  First 

Order  Model 

Mean 

Number 

Aircraft 

Coded  Values 

Response 

Casualty  Rate 

of  Acft 

Capacity 

X  Y  Z 

In  U 

1400 

14 

110 

-1 

-1 

-1 

2.894 

1400 

14 

130 

-1 

-1 

1 

2.473 

1400 

10 

110 

-1 

1 

-1 

2.435 

1400 

10 

130 

-1 

1 

1 

2.447 

1900 

14 

110 

1 

-1 

-1 

4.424 

1900 

14 

130 

1 

-1 

1 

2.854 

1900 

10 

110 

1 

1 

-1 

3.439 

1900 

10 

130 

1 

1 

1 

2.448 
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The  general  formula  for  coding  the  variables  in  this  table  is 
actually  quite  sinple.  If  L  and  U  represent  the  lower  and  upper  bounds 
of  allowable  values  for  a  particular  predictor  variable  and  if  V 
represents  the  uncoded  value,  the  coded  value  C  is  defined  as> 

C  ■  2(  V-<  L+U  )/2 )/( U-L )  or  C  ■  2<  V  -  average  )/dif ference 
For  exanple,  the  uncoded  Mean  casualty  rate  of  1600  is  equal  to  the 
coded  value  of> 

C  »  2(  1600-(  1600+1900  )/2  )/( 1900-1600  >  -  2< 1600-1750 1/300  -  -1 
The  uncoded  CRAF  capacity  of  12 6  is  equal  to  the  coded  value  oft 
C  ■  2<126-(  110+130 1/2 )/( 130-110)  -  2(  126-120  >/20  »  0.6 
Note  that  the  nid-value  of  any  uncoded  variable  is  equal  to  a  coded 
value  of  zero.  The  coded  values  Hake  the  Manual  calculations  easier  and 
enables  curve  fitting  to  be  done  using  the  Method  of  orthogonal 
polynoMials  <  25 * 229 ) j  however,  general  regression  Methods  could  also  be 
used  with  the  uncoded  variables. 


Eacaacicr.  EliiuaiifiD.  The  experinental  design  process  concluded 
that  the  effects  of  each  Main  factor,  all  two-way  interactions,  and  the 
three-way  interaction  were  significant  in  predicting  the  response 
variable  value.  A  statistical  package  ( BMDP )  was  used  to  estiwate  the 
coefficients  of  the  polynonial  terns  using  the  Method  of  least  squares 
regression.  The  least  squares  estinators  and  regression  statistics  for 
the  Model  are  provided  in  Appendix  C.  The  fitted  linear  regression 
Model  in  the  transforMed  response  variable  isi 

In  U  ■  3.033  +  0.309X  -  0.18Y  -  0.108XY 
;  -  0.326Z  -  0.261XZ  +  0.132YZ  +  0.07XYZ 

where  U  =  the  Mean  tine  in  the  systen 
X  *  coded  Mean  casualty  rate 
/  Y  *  coded  nunber  of  B -767  aircraft 

i  Z  »  coded  aircraft  capacity.  J 


The  versatility  of  this  equation  can  be  shown  in  the  following 
exanplesJ  To  find  the  naan  TIS  whan  tha  naan  casualty  rata  is  1800  par 
day  and  tha  capacity  of  aach  of  17  aircraft  is  124,  follow  tha  three- 
stap  procadurai 

1.  Coda  (or  transforn)  all  known  variablas. 

X  ■  2< 1000-1730 )/300  -  0.333 

Y  -  2(  17-17  )/2  -  0 

Z  ■  2(  124-120 )/20  -  0.4 

2.  Substituta  in  tha  aquation  and  solva  for  tha  unknown  variabla. 

In  U  ■  3.033  +  0.309<  0.333)  -  0.326(0.4)  -  0.261(0.333X0.4) 

In  U  *  2.971 

3.  Dacoda  (or  transforn)  tha  rasult. 

In  U  *  2.971 

U  -  axp( 2.971) 

U  *  19.3  hours  naan  tina  in  tha  systan 

Therefore,  undar  tha  conditions  of  tha  sinulation  nodal  and  tha 
ragrassion  techniques,  tha  patient's  avaraga  tina  in  tha  evacuation 
systan  is  19.5  hours.  Algebraic  nanipulation  of  tha  aquation's  tarns 
anablas  any  of  tha  independent  variablas  to  be  deternined  given  an 
avaraga  TIS  and  values  for  tha  other  variablas.  Thus,  if  tha  aircraft 
capacity  is  known  (say  118)  and  tha  analyst  wants  to  know  how  nany 
aircraft  will  be  required  if  tha  naan  casualty  rata  is  1680  par  day  and 
20  hours  naan  TIS  is  assunad  reasonable. 

1.  Coda  (or  transforn)  tha  variablas. 

U  -  20,  therefore  In  U  -  ln(20>  ■  2.996 

X  -  2(1630  -  1750 )/300  «  -0.467 

Z  -  2(  113  -  120  )/20  *  -0.2 

2.  Substituta  and  solva  for  Y. 

2.996  -  3.033  +  0.309( -0. 467  )  -  (0.18)Y  -  0. 108< -0.467  )Y 

-  0.326(  -0.2  )  -  0.26K -0.467  )( -0.2) 

+  0,132(  -0.2  >Y  +  0.07( -0.467  )( -0.2  )Y 

. 149Y  *  -0.067 
Y  ■  -0.449 
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3.  Decode  the  variable 


2( V  -  17  )/2  -  -0.449 

V  ■  16. SS  aircraft 

17  aircraft  will  be  required  under  the  given  conditions;  however, 
this  example  deHonstrates  the  actual  discontinuity  of  the  independent 
variables  even  though  the  response  surface  is  treated  as  a  continuous 
surface.  This  does  not  cause  problems  as  long  as  the  results  are 
interpreted  correctly. 

Because  the  number  of  aircraft  and  the  aircraft  capacity  can  assume 
only  integer  values,  the  response  surface  results  must  be  rounded  when 
solving  for  these  variables. 

QlbftE  Bcspauac  Surfaces-  Three  additional  response  surface 
equations  were  generated  using  the  same  regression  techniques;  however 
additional  simulation  experiments  were  required.  The  regression  results 
are  provided  in  Appendix  C  for  the  response  surfaces  shown  in  Table 
VIII .  The  same  procedures,  outlined  previously,  apply;  but,  one  must 
now  have  some  idea  of  the  number  of  aircraft  operating  in  the  system. 

Table  VIII 

Linear  Response  Surfaces 

Response  Surface  Equation 
In  U  ■  3.212  +  0.429X  -  0.4482  -  0.337XZ 

In  U  •  3.028  +  0.336X  -  0.343Z  -  0.303XZ 

In  W  -  2.861  +  0.202X  -  0.201Z  -  0.191XZ 

i 

Using  the  same  data  as  in  the  first  example  <1800,17  and  124),  the 
response  equation  corresponding  to  17  aircraft  yields  a  mean  TIS  ofs 
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16 

17 

13 
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In  U  ■  3.028  ♦  0.334(  0.333 )  -  0.343(0.4)  -  0.303(0.333X0.4) 

In  U  *  2.942,  or  II  *  19.3  hours 

The  previous  exanple  resulted  in  «  rte«n  T1S  of  19.5  hours.  To  see 
how  both  the  two-  end  three-verieble  response  surface  equetion  conpare 
with  sinulation  experihents,  assune  the  following  paraneter  velues  for 
sinplicityi 

ween  cesuelty  ret#  ■  1750  or  coded  velue  *  0 

nunber  of  eircreft  ■  14  or  coded  velue  ■  -1 

eircreft  cepecity  ■  130  or  coded  velue  »  +1 

Three-verieble  Model  yielded  e  Mean  TIS  of  15.7  hours  end  the  two- 
verieble  Model  yielded  e  Mean  TIS  of  15.9  hours.  A  three  replicete 
sinuletion  everage  yielded  e  Mean  TIS  of  14.8  hours.  An  interpretation 
in  ter  ms  of  the  evacuation  systeM  will  be  provided  in  the  next  chapter; 
however,  the  sonewhet  higher  predicted  velues  indicate  that  a  linear 
exponential  expression  May  not  be  an  appropriate  functional  forM  of  the 
response  surface. 

Quadratic  Bedels 

In  order  to  develop  a  quadratic  functional  forM  for  the  exponential 
Model,  sinuletion  experiMents  were  conducted  at  the  Midpoint  values  of 
each  factor.  Now  however,  since  the  design  is  no  longer  orthogonal,  the 
coding  of  high  and  low  values  is  not  appropriate. 

The  only  coding  perforned  in  developing  a  quadratic  tiodel  was  to 
divide  the  casualty  rate  factor  by  100  and  the  aircraft  capacity  by  10. 
An  exaMple  of  a  typical  data  point  with  this  coding  night  be  (14,18,11) 
corresponding  to  a  casualty  rate  of  1400  and  18  aircraft  with  a  capacity 
of  1?0. 


Eicioaitc  At  in  the  linear  model,  a  BMDP  statistical 


package  was  used  to  estimate  the  coefficients  of  the  quadratic 
polynomial.  An  "all  subsets"  regression  was  performed  using  all  linear 
and  quadratics  wain  effects  and  all  two-way  interactions.  The  least 
squares  estimators  and  regression  statistics  for  the  selected  "best" 
subset  are  shown  in  Appendix  C.  The  following  quadratic  regression 
model  was  determined! 

In  M  -  -2.98  +  2.095X  -  1.81Z  +  0.0613XX  +  0.0966ZZ 

-  0.1078XY  -  0.1833XZ  +  0.1408YZ 

where  U  ■  mean  time  in  the  system 
X  ■  casualty  rate/100 
Y  »  number  of  aircraft 
Z  *  aircraft  capacity/10. 

The  quadratic  model  has  the  same  versatility  as  the  linear  model, 
but  yields  response  surface  values  closer  to  the  experimental  results. 
Th*  same  treatment  combinations  as  in  the  linear  model  will  be  used  to 
compare  the  models.  If  the  casualty  rate  is  1800  per  day  and  each  of  17 
aircraft  has  a  capacity  of  124,  the  mean  time  in  the  system  (TIS)  is 
found  by  encoding,  substituting  and  solving  for  U. 

In  U  -  -2.98  ♦  2.093(18)  -  1.81(12.4)  +  0.0613(  18 )( 18  > 

+  0.0966(  12.4X12.4)  -  0.1078(  18X17) 

-  0.1833(18X12.4)  +  0.1408(17X12.4) 

In  U  -  2.7317 

U  *  16.15  hours  mean  TIS 

QibftC  SespfiQse  Surfaces •  Again,  as  in  the  linear  case,  additional 
response  surface  equations  were  developed  for  the  cases  where  the  number 
of  aircraft  is  known.  Appendic  C  also  contains  the  least  squares 
estimators  and  regression  statistics  for  these  surfaces.  Table  IX 
summarizes  the  response  surface  equations. 
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Table  IX 


Quadratic  Response  Surfaces 


Aircraft 

Response  Surface  Equations 

14 

In  U  - 

3.100  *  0.0810XX  +  0.13427.Z  -  0.2108XZ 

17 

In  U  - 

2.751  +  0.0737XX  ♦  0.1316ZZ  -  0.1978XZ 

18 

In  U  - 

2.544  +  0.0441XX  +  O.O840ZZ  -  0.1248XZ 

Model  Cfiopaciifica 

For  the  particular  treatnent  conbination  of  casualty  rate  -  1750, 
and  each  of  17  aircraft  with  capacity  of  110,  the  five  nodels  (including 
a  three-replicate  average  of  the  sinulation  Model )  produced  the 
following  Mean  TIS  values* 


Model 

Mean 

TIS 

SiHulation  ExperiMent 

22.8 

hours 

Three-variable  linear  Model 

28.8 

hours 

Two-variable  linear  Model 

29.1 

hours 

Three-variable  quadratic  Model 

24.7 

hours 

Two-variable  quadratic  Model 

23.8 

hours 

The  quadratic  Models  are  clearly  More  accurate  than  the  linear 
Models.  Interpretation  of  these  results  in  terMS  of  the  evacuation 
systeH  will  be  included  in  chapter  V.  The  functional  forn  of  the 
response  surface  equation  that  More  accurately  reflects  the  experiMental 
results  has  been  established.  The  results  of  the  exponential  equation 
with  a  quadratic  polynoxial  exponent  can  be  used  to  provide  insight  into 
the  questions  surrounding  an  "acceptable"  evacuation  and  distribution 
systeM. 
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Applying  response  surface  nethodologies  to  the  casualty  avacuation 
sinulation  nodal  has  providad  an  afficiant  Method  for  ansuaring  nany 
questions  partaining  to  tha  affacts  tha  thraa  indapandant  variables  have 
in  tha  nodal.  Tha  exanples  prasantad  in  this  chaptar  ara  sona  of  tha 
ways  in  which  a  rasponsa  surfaca  can  ba  usad  to  halp  analyza  tha  real- 
world  casualty  avacuation  problan. 

Tha  next  chaptar  will  examine  tha  rasults  of  tha  sinulation 

axparinants  and  tha  rasponsa  surfacas  in  tarns  of  tha  systan  it 

sinulatas  or  dascribas.  Tha  varsatility  of  tha  sinulation  nodal  will  ba 

danonstrated  and  sona  obsarvations  on  tha  systan' s  structura  and 

parfornanca  will  ba  providad.  Conclusions  and  raconnandations  will 
« 

follow  in  Chaptar  VI. 
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V.  EXPERIMENTAL  RESULTS 

loi.CQduci.ifio 

i 

The  research  effort  hat  progressed  froH  the  developMent  of  a 
casualty  avacuation  siMulation  Modal ,  to  an  axparinantal  design  procass 
and  finally  to  ratponsa  surface  Methodology.  Each  step  of  the  sequence 
has  produced  various  kinds  of  data  and  relationships.  This  chapter  will 
provide  observations  and  explanations  of  the  results  fron  these 
developMents. 

The  basic  evacuation  structure  has  been  outlined  in  Chapter  II  and 
included  the  range  of  the  independent  variables  and  specific  values  for 

the  siMulation  experinents.  Chapter  III  established  the  statistically 

j 

significant  variables  and  fixed  the  values  for  other  controllable 
factors.  Response  surfaces  were  generated  in  Chapter  IV. 

The  analysis  of  results  includes  the  siMulation  experinents 
conducted  for  sone  pilot  runs  and  for  those  experinents  that  provided 
the  response  values  used  in  the  experinental  design  and  response 

i 

surfaces.  Other  explanations  will  include  results  fron  variations  Made 
to  a  base  case  Model. 

Siuulatifio^E^peciaefiisl  Design  Besulis 

In  addition  to  providing  a  Means  of  Model  verification,  prelininary 
sinulation  runs  also  helped  establish  a  feeling  for  various  Measures  of 
evacuation  perfornance  and  a  range  of  values  for  then.  One  early 
experinent  did  not  differentiate  between  those  patients  offloading  at 
the  staging  bases,  those  continuing  to  other  U.S.  cities,  and  those 
patients  causing  an  overflow  situation. 
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Though  th*  overall  titan  tint  in  tha  system  (TIS)  remained  the 
primary  measure  of  avacuaiion  performance,  tha  number  of  casualties  in 
aach  catagory  abova  and  tha  maximum  TIS  for  aach  catagory  wara  important 
aacondary  measures  that  cauiad  aona  nodal  nodificationa.  Originally, 
all  aircraft  wara  required  to  ba  fully  loadad  bafora  dapartura.  In  some 
caaaa,  thia  caused  an  axcaaaiva  anount  of  overflows  therefore,  C-9 
aircraft  are  now  permitted  to  depart  at  752  capacity  if  lass  than  a  full 
load  of  casualtias  are  waiting  for  transportation. 

The  pilot  runs  also  pointed  out  specific  treatment  combinations 
resulting  in  grossly  unacceptable  mean  TIS  values.  These  results 
established  bounds  for  the  factors  used  in  the  experimental  design 
process. 

B~Z6Z  Utilization.  Though  the  study  did  not  intend  to  maximize 
B-767  utilization  or  optimize  the  number  of  aircraft  required,  high 
utilization  rates  were  achieved  in  all  but  one  treatment  combination. 

A  combination  of  low  casualty  rate  and  high  capacity  with  18  operating 
aircraft  yielded  a  7 OX  utilization  rate.  All  other  design  combinations 
resulted  in  utilization  rates  greater  than  802.  Overcommittment  of 
resources  occurred  when  the  system  was  saturated.  This  happened  at  only 
two  design  points  -  a  high  casualty  rate  and  low  capacity  combination 
with  either  16  or  18  aircraft. 

C_2  UiiliZalifiD.  Though  the  experimental  design  process  concluded 
that  the  number  of  C-9  aircraft  was  not  statistically  significant,  the 
average  C-9  utilization  rate  was  between  70X  and  802  when  only  six  C-9 
aircraft  were  exercised.  When  nine  C-9  aircraft  were  used,  the  rate 
ranged  from  502  to  652.  Six  C-9  aircraft  might  be  sufficient  to 
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I  distribute  casualties  if  an  efficient  CRAF  scheduling  policy  is  adopted 

and  if  potential  patient  overflow  is  not  increased. 


| 


Such  a  policy  would  schedule  a  CRAF  aircraft  arrival  only  if  at 
least  three  C-9  aircraft  were  available.  This  would  maximize  C-9 
aircraft  and  aircrew  utilization  and  minimize  patient  overflow 
situations  described  in  chapter  II.  Deternining  a  "best"  scheduling 
policy  as  well  as  the  optimum  number  of  C-9  aircraft  would  require  an 
additional  study. 

ttiKitma  list  is  ib*  Syaiea 

The  Meaning  of  maximum  TIS  should  be  explained  again  for  emphasis. 
It  includes  the  tine  a  patient  spends  in  flight,  the  ground  tine  between 
aircraft  transfers  when  necessary,  the  ground  tine  for  CRAF  pre¬ 
departure  activities  and  the  tine  spent  in  the  hospital  awaiting 
transportation. 

The  actual  maximum  amount  of  tine  spent  outside  of  the  formal 
hospital  environment  for  any  patient  is  less  than  18  hours.  A  two-hour 
bus  or  anbulance  ride  at  final  destination  would  increase  this  tine  to 
20  hours.  The  expected  tine  in  transit  is  about  14  hours)  therefore, 
the  maximum  amount  of  tine  spent  waiting  for  transportation  is  14  hours 
less  than  the  tiaxiHUH  TIS. 

bean  and  battiaud  IIS 


ft 


Except  in  those  situations  with  very  little  patient  queueing,  the 
maximum  TIS  value  is  approximately  twice  the  mean  TIS  value.  This  is 
partly  due  to  the  first-in/first-out  (FIFO)  queue  discipline  of  the 
model.  Also,  evacuation  occurs  at  a  fairly  constant  rate)  therefore, 
the  max  TIS  should  be  about  twice  the  the  mean  TIS. 
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The  convenient  relationship  between  the  Mean  and  naxinuM  TIS 
provides  a  useful  Method  for  deternining  an  acceptable  value  for  the 
priMary  response  variable.  This  will  be  conbined  with  the  other 
Measures  of  systeH  perforxance  to  establish  a  realistic  level  of  the 
Mean  TIS. 
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CdiecifiD  Value  o£  btso  IIS 

Selecting  a  value  for  the  Mean  TIS  will  depend  on  its  inplications 
on  whether  the  evacuation  structure  can  handle  the  corresponding 
workload  and  if  the  values  of  the  secondary  Measures  of  Merit  "wake 
sense".  Obtaining  realistic  values  would  require  additional  studies; 
however,  for  this  research,  subjective  values  will  be  used  to  sinplfy 
the  probleH  and  to  show  how  these  secondary  Measures  can  be  used  to  help 
establish  a  useful  criterion  value  of  Mean  TIS. 

In  addition  to  the  Mean  TIS  value,  other  inportant  Measures  of 
perforMance  include  the  MaxixuM  TIS  and  the  MaxinuM  nunber  of  casualties 
awaiting  airlift.  Mean  waiting  tine,  nunber  of  overflow  patients  and 
average  nunber  of  waiting  patients  Must  also  be  considered.  Each 
analyst  Must  deternine  a  Mean  TIS  criterion  value  to  judge  treatnent 
conbinations.  What  follows  is  one  possible  developnent  an  analyst  can 
pursue  and  indicates  sone  factors  that  Must  be  evaluated  in  selecting  a 
criterion  value. 

Based  on  the  initial  casualty  distribution  assunption  that  the 
highest  casualty  rate  is  2S00  per  day,  a  naxinun  of  2S00  patients 
waiting  for  transportation  will  be  considered  acceptable.  If  actual 
evacuation  operations  begin  no  later  than  24  hours  after  the  conflict 
starts,  the  five  echelons  of  care  should  be  capable  of  absorbing  this 
Many  casualties. 
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A  reasonable  naxinun  waiting  tin*  for  a  patiant  could  be  36  hours. 
Whan  tha  expected  in-transit  tine  of  14  hours  is  addad  to  this  value,  a 
naxinun  TIS  of  SO  hours  would  not  be  unraalistic.  Basad  on  tha  ganaral 
.  *  ralationship  batwaan  naan  and  naxinun  TIS,  a  starting  point  for 

,  accapting  or  rajacting  spacific  treatnent  conbinations  would  ba  a  naan 

TIS  valua  of  25  hours. 

If  tha  nunbar  of  ovarflow  patiants  listad  in  Tabla  I  was  not 
axcaadad  and  tha  othar  naasuras  of  parfornanca  ara  accaptabla,  a 
rasponsa  valua  of  25  hours  would  inply  an  accaptabla  factor  conbination. 
By  conbining  all  accaptabla  factor  conbinations,  dacisions  concarning 
tha  nunbar  of  aircraft  to  contract  or  how  larga  a  naan  casualty  rata  tha 
systan  could  support  can  ba  nada.  A  tabla  of  thasa  results  will  ba 
provided  later. 

Using  tha  25-hour  criterion  causes  two  of  the  axperinental  design 
treatnents  to  ba  elininated.  A  high  casualty  rata  and  low  aircraft 
capacity  conbination  results  in  larga  naan  TIS  values,  excessive  naxinun 
TIS  values,  large  nunbers  of  waiting  casualties  and  a  large  overflow 
population.  Many  other  experinental  runs  ara  required  to  deternina  tha 
naan  TIS  for  values  between  tha  low  and  high  linits  of  each  factor  in 
tha  design.  Tha  valua  of  tha  response  surfaces  is  clearly  shown  herej 
*  •  however,  sone  interior  values  can  still  be  sinulated. 

flidcaiDi  Eapcriacuis 

This  section  describes  the  results  of  the  experinents  conducted  at 
|  nidpoint  values  of  the  independent  variables.  Specifically,  results 

fron  conbinations  of  17  CRAF  aircraft  with  a  120-patient  capacity  and  a 
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casualty  rata  of  1750  par  day  will  ba  examined  for  data  consistancy  and 
conparison  with  tha  linaar  rasponaa  surfaca  pradictions. 

Tha  data  in  Tabla  X  summarizes  all  tha  treatment  combinations 
possibla  using  at  least  ona  Midpoint  value.  Each  antry  is  an  avaraga  of 
thraa  raplicatas  of  tha  simulation  nodal.  Tha  raw  data  was  usad  to 
davalop  tha  quadratic  rasponsa  surfaca  aquations. 

Tabla  X 


Sample  Maan  TIS  (  Hours  > 


Casualty  Rate 

1400 

1750 

1900 

B-767 

Capacity 

110 

120 

130 

110 

120 

130 

110 

120 

130 

16  CRAF 

18.2 

15.0 

14.5 

40.9 

17.5 

14.8 

83.7 

42.4 

17.5 

17  CRAF 

15.4 

14.6 

14.2 

22.8 

15.3 

14.6 

55.3 

24.6 

15.3 

18  CRAF 

14.4 

14.1 

13.9 

15.5 

14.6 

14.2 

31.4 

17.6 

14.4 

As  expected,  tha  increasing  or  decreasing  trends  were  verified  with 
these  intermediate  data  points.  This  relationship  can  be  useful  in 
selecting  or  eliminating  certain  treatment  combination  easily.  For 
example,  since  the  average  of  three  simulation  experiments  with  the  same 
treatment  combination,  say  ( 1750, 17, 110  >,  is  less  than  25  hours,  then 
tha  result  of  decreasing  tha  casualty  rate,  increasing  the  number  of 
aircraft,  or  increasing  the  capacity  would  be  less  than  25  hours.  If 
all  other  measures  of  performance  were  satisfactory,  these  "better" 
treatment  combinations  would  also  be  acceptable. 

It  is  useful  to  hold  two  of  the  independent  variables  fixed  and 
observe  the  rate  of  change  of  tha  response  variable  as  the  third 
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independent  variable  change*.  In  som#  cases,  the  Mean  TIS  appears  to 
have  a  linear  relationship,  but  a  very  non-linear  relationship  in  other 
cases.  This  intuitive  relationship  strengthens  the  developnent  of 
quadratic  response  surfaces. 

Casualty  Bata  Ebpariaasts 

The  triangular  distribution  assuned  at  the  beginning  of  the  study 
was  changed  to  see  how  other  distributions  affect  the  Measures  of 
perfornance.  The  base  case  triangular  distribution  had  parameter  values 
of  <1000,1750,2500)  resulting  in  a  Mean  TIS  of  16.7  hours  in  a  systeM 
with  16  aircraft  with  capacity  of  120.  A  Unifornt 1000,2500)  and  a 
NornaK  1750,250)  were  tested  under  the  sane  conditions  resulting  in  19.9 
hours  and  15.6  hours  respectively. 

The  apparent  trend  indicates  that  a  distribution  with  a  larger 
standard  deviation  results  in  a  larger  Mean  TIS  when  all  other  factors 
are  held  constant.  In  these  cases,  the  standard  deviations  range  fron 
250  for  the  Norrtal  distribution  to  433  for  the  Unifornf  the  triangular 
distribution  standard  deviation  is  306. 

Another  variation  in  the  casualty  rat*  consisted  of  a  cyclic 
generation  of  casualties.  This  Might  approxinate  the  cycle  between 
heavy  conbat  losses  and  regeneration  periods  with  few  casualties.  A  10- 
day  cycle  between  successive  high  casualty  rates  beginning  at  a  Mean 
value  of  1750  per  day  was  exaMined.  The  sinusoidal  function  used  to 
siMulate  the  cycle  is  listed  here  and  in  the  conputer  code. 

Casualty  rate  *  1750  +  750*SIN< 2*PI*D0U/10 ) 
where  DQU  represents  the  day  of  the  war. 
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This  variation  resulted  in  a  naan  TIS  of  24.6  hour*,  an  increase  of 

i 

eight  hour*.  The  patient  queueing  that  rasults  fron  haavy  casualty 

I 

rata*  raquirad  about  two  da^s  of  intansiva  evacuation  operations  to 
relieve  the  hospital  congestion.  The  FIFO  queue  discipline  causes  a 
greater  waiting  tine  for  those  casualties  generated  on  the  "down-hill" 
portion  of  the  cycle.  However,  the  cyclic  generation  of  casualties 
would  probably  be  danpened  sonewhat  as  the  casualties  process  through 
the  four  or  five  echelons  of  care  before  evacuation. 

Ibc  22-Uque  Criterion 

The  level  of  this  criterion  was  established  subjectively,  but  based 
on  reasonable  levels  of  other  neasures  of  perfornance.  Sone  treatnent 
conbinations  that  would  be  rejected  according  to  this  criterion  night 
still  considered  acceptable  in  certain  instances.  This  gray  area 
consists  of  response  values  slightly  greater  than  25  hours,  by  one  or 
two  hours  for  instance. 

One  particular  case  had  a  nean  TIS  of  26.5  hours  but  at  the  30-day 
point  in  the  war,  the  nean  TIS  was  25  hours  and  all  other  Measures  of 
nerit  were  within  acceptable  linits.  A  look  at  the  systen  at  the  30-day 
point  is  well  after  the  initial  surge  period  and  sea  lines  of  resupply 
will  have  been  established.  A  decreasing  C-141  resupply  requirenent 
suggests  that  if  additional  airlift  could  be  connitted  at  that  tine, 
additional  treatnent  conbinations  could  becone  feasible  alternatives. 


EatieDi  Becsyery  Iitie  / 

The  base  case  nodel  <1750,16,120)  included  a  specific  anount  of 
hospital  recovery  tine  for  each  category  of  patient.  Another  variation 


of  this  nodel  included  increasing  the  recovery  periods  by  25%. 


The 
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result  was  an  increase  of  the  naan  TIS  to  20.1  hours.  This  reflects 
those  situations  when  a  hospital  is  full  and  remains  full  for  a  longer 
period  thus  requiring  More  and  longer  flights  to  other  hospitals. 

Minor  overflow  conditions  were  recorded  at  two  of  the  staging  bases 
that  previously  had  none.  The  increased  recovery  tines  adversly 
affected  the  maximum  TIS  and  maximum  number  of  waiting  casualties  to  the 
point  where  the  treatment  conbination  would  be  rejected.  The  maximum 
TIS  increased  from  31.3  to  56.3  and  the  maximum  number  in  the  waiting 
queue  increased  from  1080  to  2760.  Though  there  is  no  reason  to  suspect 
the  recovery  tine  estimates,  a  marginally  acceptable  option  could  cause 
serious  problems  if  the  estimates  are  too  low. 

Bcspqdsa  Surface  BcsuUi 

To  avoid  the  costly  and  time  consuming  process  of  making  additional 
simulation  experiments,  response  surface  equations  were  developed  in 
chapter  IV.  Since  the  experimental  design  data  points  were  used  to 
determine  the  coefficients  of  the  terms  in  the  linear  equations,  the 
simulation  results  should  be  closely  predicted  by  the  linear  response 
surface  equations  evaluated  at  those  design  points. 

Simulation  experiments  were  conducted  at  all  combinations  of 
midpoint  values;  quadratic  response  surface  equations  were  developed 
from  the  results.  The  data  in  Table  XI  summarizes  the  average 
simulation  results  and  predicted  values  from  the  three-variable 
quadratic  equation.  The  error  term  measures  the  deviation  from  the 
predicted  value. 
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Table  XI 


Comparison  of  Simulation  and  Response  Surface  Results 


Treatnent 

Predicted 

Sinulated 

Error 

1600,16,110 

18.8 

18.2 

-  32 

,120 

16.3 

15.0 

+  52 

,130 

14.8 

15.3 

♦  9% 

1600,17,110 

15.7 

15.4 

-  22 

,120 

13.8 

14.6 

+  62 

,130 

14.8 

14.2 

-  42 

1600,18,110 

13.2 

14.4 

*  92 

,120 

13.4 

14.1 

♦  52 

,130 

16.4 

13.9 

-152  * 

1750,16,110 

34.6 

40.9 

•*■182  * 

,120 

20.1 

17.5 

-132  * 

,130 

14.1 

14.8 

+  52 

1750,17,110 

24.7 

22.8 

-  82 

,120 

16.5 

15.3 

-  72 

,130 

13.4 

14.6 

+  92 

1750,18,110 

17.6 

15.5 

-122  * 

,120 

13.5 

14.6 

♦  82 

,130 

12.6 

14.2 

♦132  * 

1900,16,110 

84.1 

83.7 

-  12 

,120 

37.1 

42.4 

+142  * 

,130 

19.8 

17.5 

-122  * 

1900,17,110 

51.0 

55.3 

+  82 

,120 

25.9 

24.6 

-  52 

,130 

15.9 

15.3 

-  42 

1900,18,110 

31.0 

31.4 

+  12 

,120 

18.1 

17.6 

-  32 

,130 

12.8 

14.4 

+132  * 

Sone  of  the  errors  nay  s.een  unacceptably  large  in  the  table,  but 
note  the  regions  where  the  larger  errors  occur  in  terns  of  the  nean  TIS 
value.  The  largest  errors  correspond  to  treatnents  whose  predicted  nean 
TIS  is  below  20  hours  or  above  30  hours  approxinately.  Based  on  the 
criterion  value  deternined  earlier  (25  hours),  the  decision  regarding 


whether  a  treatment  combination  ia  accepted  or  rejected  remains 
unchanged.  Thus,  the  (1750,16, 110 )  conbination  is  rejected  even  though 
the  prediction  underestinates  the  sinulated  value}  the  ( 1900, 16, 130 ) 
combination  is  accepted  even  though  the  prediction  overestimates  the 
simulated  value. 

If  an  analyst  selects  a  criterion  value  below  20  hours  or  above  30 
hours,  the  problem  would  require  re-examination  of  the  evacuation 
structure  and  parameter  values.  Based  on  this  study's  structure  and  its 
measures  of  performance,  such  high  or  low  criterion  values  appear 
unlikely. 

Treatment  combinations  with  predicted  mean  TIS  values  near  25  hours 
are  within  102  of  the  experimental  averages.  Is  this  degree  of 
predictive  accuracy  close  enough?  In  an  advanced  planning  situation 
with  many  assumptions  and  unknowns,  the  quadratic  model  can  be  used  with 
confidence.  It  adequately  measures  the  performance  of  a  system 
described  by  the  casualty  rate,  number  of  aircraft  and  aircraft 
capacity  and  operating  according  to  this  study's  structure. 

The  response  surface  equations  generally  predict  larger  values  than 
this  particular  sample.  In  other  tests,  the  same  relationship 
prevailed.  The  built-in  conservative  nature  of  the  response  surface 
equations  provides  a  margin  of  safety  for  selecting  viable  treatment 
combinations. 

EiDal  Besulia 

The  examples  in  Chapter  IV  provide  some  idea  of  the  kinds  of 
questions  that  can  bo  answered  by  using  the  response  surface  equations. 
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Given  the  Model's  parameters  in  this  study  end  the  23-hour  criterion 
value,  the  data  that  follows  in  Table  XII  will  help  answer  the  likely 
questioni  Uhat  are  the  upper  limits  of  the  Mean  casualty  rate  that  a 
systeM  coMprised  of  X  aircraft  with  capacity  Y  can  evacuate? 


Table  XII 

Three-variable  Quadratic  Response  Surface  SuMHary 


Number  of  B-767 

B-767  capacity 

Casualty  rate 

16 

110 

1675 

16 

112 

1700 

16 

114 

1725 

16 

116 

1755 

16 

118 

1780 

16 

120 

1810 

16 

122 

1840 

16 

124 

1870 

16 

125 

1885 

16 

up  to  130 

1900 

17 

no 

1750 

17 

112 

1775 

17 

114 

1805 

17 

116 

1830 

17 

118 

1860 

17 

120 

1890 

17 

up  to  130 

1900 

18 

110 

1850 

18 

112 

1875 

18 

113 

1890 

18 

up  to  130 

1900 

Values  greater  than  a  Mean  of  1900  casualties  per  day  could  be 
evacuated  but  additional  studies  should  be  Made  on  the  tactical 
evacuation  capabilities  and  the  stateside  distribution  requirements. 
The  table  shows  that  for  a  fixed  number  of  aircraft,  the  Mean  TIS 
increases  by  about  15  casualties  per  day  per  unit  increase  in  aircraft 
capacity.  Each  aircraft  is  carrying  one  extra  patient  each  day, 
approximately. 


Sutmacy 


The  siMulation  experiMents  provided  various  Measures  of  systen 
perforMance.  High  aircraft  utilization  rates  were  achieved  at  Most 
operation  conditions.  With  an  intelligent  scheduling  algorithM  for  the 
CRAP  operations,  a  C-9  operation  with  as  few  as  six  aircraft  May  be 
achievable. 

A  realistic  MaxiMUM  nuxber  of  waiting  casualties  was  established 
and  a  general  relationship  between  Mean  and  MaxiMUM  TIS  was  observed. 
These  provided  an  initial  estiMate  of  a  decision  or  criterion  value  of 
the  Mean  TIS. 

This  criterion  value  established  a  basis  for  coMparing  other 
treatMent  coMbinations.  The  effects  of  varying  original  paraMeters  of 
the  Model  (values  and  distributions)  were  judged  according  to  the 
criterion  value. 

SiMulation  results  were  coMpared  to  the  quadratic  response  surface 
predictions.  The  quadratic  equations  predicted  values  closer  to  the 
experiMental  results  than  the  linear  response  surface  equations. 

The  suMMary  of  results  based  on  the  Model's  assuMptions,  correction 
factor  and  criterion  value  for  the  Mean  TIS  provides  a  guide  for 
evaluating  the  perforMance  and  structure  of  the  evacuation  systeM 
envisioned.  Conclusions  and  recoMMendations  will  follow  in  the  last 


chapter 


VI.  CONCLUSION 


Btwitw 

The  Casualty  Evacuation  Modal  ( CEM )  developed  in  this  research 
affort  began  with  tha  realization  that  tha  outcona  of  a  conventional 
Europaan  conflict  could  ba  affectad  by  tha  corniitnent  of  airlift 
rasourcas  to  tha  avacuation  Mission.  Currant  studies  ara  evaluating  tha 
usa  of  Hodifiad  B -767  aircraft  for  casualty  avacuation.  This  study  has 
attenpted  to  answsr  tha  question  of  how  an  avacuation  and  distribution 
natwork  using  B -767  and  C-9  aircraft  can  ba  structurad  to  parforw  tha 
task  and  how  its  perfornance  can  ba  avaluatad.  A  sinulation  nodal  was 
davalopad  tc  provide  sons  answers  and  insight  to  tha  problan. 

SlEUCiUCft 

Though  a  60-day  conflict  was  used  for  tha  study,  tha  structure  of 
tha  natwork  is  essentially  independent  of  tha  conflict's  duration.  Tha 
structura  depends  on  casualty  rates,  nunbers  and  types  of  aircraft 
involved,  tha  nunber  and  location  of  oparating  basas,  and  Many  other 
factors. 

Various  casualty  rates  were  examined  in  the  study  ranging  fron  1000 
to  2300  par  day.  Two  types  of  aircraft  were  considered  in  the 
avacuation  plan,  while  the  operating  bases  included  one  onload  location, 
four  staging  bases  and  wore  than  40  offload  locations  in  the  eastern 
half  of  the  U.S. 

Tha  B-767  aircraft  performed  the  irierlheater  evacuation  task  and, 
when  necessary,  provided  some  distribution  services.  The  C-9  aircraft 
operated  frow  one  location  interfacing  with  the  B-767s.  With  these 


airlift  resources,  MAC  would  provide  the  link  between  the  BoD  Theater 
CoMbat  Medical  System  and  the  National  Disaster  Medical  System. 

Four  staging  bases  on  the  eastcoast  -  Boston,  New  York, 
Philadelphia  and  Washington  -  provide  access  to  More  than  25,000 
hospital  beds  in  More  than  25  cities  that  May  participate  in  the  NDMS. 
The  C-9  network,  operating  froM  Washington,  provides  access  to  More  than 
35,000  additional  beds  in  30  cities. 

CftEfrttBBQCft 

The  priMary  Measure  of  the  system's  perforMance  is  the  aMount  of 
tine  a  patient  spends  in  the  systeM.  This  tine  includes  all  flight 
tiMes,  predeparture  ground  tiMes,  and  the  tiMe  spent  in  the  hospital 
waiting  for  transportation.  Secondary  Measures  of  perforMance  include 
aircraft  utilization,  average  and  HaxiMUM  waiting  tiMes,  average  and 
HaxiMUM  nuxber  waiting,  and  the  nuxber  of  overflow  patients.  Overflow 
patients  are  described  in  chapter  11. 

A  value  of  25  hours  time  in  the  systeM  (TIS)  was  selected  as  the 
criterion  value  between  accepting  or  rejecting  a  systeM's  perforMance. 
The  systeM  is  identified  by  its  treatMents  ( casualty  rate,  nuMber  of 
aircraft  and  aircraft  capacity  operating  within  the  structure  outlined 
above ). 

FroM  the  results  of  selected  specific  treatment  combinations  in  the 
siMulation  model,  response  surfaces  were  developed  to  predict  a 
perforMance  value  for  treatxent  combinations  that  could  not  be 
simulated.  The  equations  developed  were  conservative  predictors  of  the 
mean  TIS  and  thus  provide  a  built-in  Margin  of  safety  in  estimating  a 
system's  performance. 


Esctfici 


The  simulation  model  utti  six  factor*  to  represent  a  tpacific 
•vacuation  network.  Tha  primary  factors  included  tha  casualty  rata, 
number  of  6 -767  aircraft  and  tha  aircraft  capacity.  Pradapartura  ground 
tines,  flight  tines  and  tha  nunbar  of  C-9  aircraft  make  up  tha  other 
factors.  All  of  tha  factors  can  ba  modified  to  accommodate  many 
variations  within  tha  structure.  Specific  levels  used  in  this  study  for 
each  factor  are  described  in  chapter  III. 

Suoascy  a£  Beaulia 

Tha  research  objectives  in  chapter  I  have  bean  achieved.  A  routing 
structure  was  designed  and  tasted  in  a  simulation  modal  under  various 
operating  conditions.  Tha  structure  included  the  number  and  location  of 
the  necessary  intermediate  staging  bases.  For  each  factor  combination, 
the  number  of  B -767  aircraft  required  can  be  estimated  and  other 
measures  of  performance,  particularly  the  mean  time  each  patient  spends 
in  the  system,  can  be  determined. 

Two  types  of  model  variations  were  examined.  The  first  variations 
allowed  changes  to  the  three  factors  affecting  the  primary  measure  of 
performance.  The  second  type  of  variation  examined  the  effects  of 
changes  made  to  the  model's  assumptions  and  raw  data  estimates.  The 
results  from  all  variations  were  compared  with  a  selected  criterion 
value  of  25  hours  mean  time  in  the  system. 

Specific  values  for  each  factor  combination  and  variation  are  not 
as  important  as  the  degree  and  direction  of  change  in  the  response 
variable.  Generally,  at  low  casualty  rates,  none  of  the  variations 
caused  an  appreciable  change  in  the  mean  TIS.  At  high  casualty  rates 


and  corresponding  high  Mean  TIS  values,  variations  did  not  inprove  the 
value  of  Mean  TIS  enough  to  cause  a  previously  rejected  factor 
conbination  to  be  accepted. 

Factor  coMbinations  and  variations  resulting  in  a  Mean  TIS  between 
20  and  30  hours  constitute  the  gray  area)  however,  the  results  have  a 
consistent  trend.  A  sinulation  experinent  yielding  25  hours  or  less  was 
found  to  be  acceptable  in  all  other  Measures  of  perfornance  as  well.  A 
response  surface  prediction  of  30  hours  or  less  also  had  acceptable 
perfornance  values. 

BtCQtmendatiQDS 

The  recoHMendations  for  further  study  deal  prinarily  with  the 
structure  of  the  evacuation  and  distribution  network.  This  study 
exaMined  one  feasible  structure  out  of  Many  possibilities  and  within 
this  structure,  other  variations  should  be  considered. 

Excluding  the  tactical  evacuation  requirenents,  the  first  area  that 
warrants  additional  attention  is  whether  Multiple  onload  locations 
should  be  incorporated  into  the  siMulation  Model.  While  flying  tiMe 
could  be  incorporated  easily,  specifying  the  nunber  and  location  of 
additional  onload  bases  and  estinating  the  nunber  of  casualties 
departing  each  location  could  introduce  More  error  than  a  single  onload 
assunption.  The  Multiple  onload  forMUlation  also  suggests  an  atteMp  to 
optiMize  the  routing  structure. 

The  choice  of  B-767/C-9  interface  location  should  also  be 
investigated.  Washington  was  selected  because,  of  the  four  staging 
bases,  it  was  closer  to  the  geographical  center  of  the  cities  served  by 
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the  C-9  fleet.  This  city  however,  had  a  greater  potential  for  a 
patiant  ovarflow  problem 

Thasa  reconnendations  suggest  fornulating  a  scheduling  policy  to 
provida  intalligant  intarfaca  batwaan  tha  two  typas  of  aircraft. 
Without  a  raalistic  plan,  tne  potantial  for  ovarflow  would  incraasa. 

Tha  study  suggastad  a  102  upward  adjustnent  to  tha  nunber  of 
aircraft  raquirad  to  conpensate  for  tha  lack  of  naintenance  and  waather 
dalays  and  divarsions.  If  thasa  factors  could  instaad  ba  addad  to  tha 
siiiulation  nodal,  it  would  nore  closaly  raprasant  tha  sinulated  real- 
world  oparation. 

Eioal  Bcback 

An  optinal  evacuation  plan  would  ba  good  to  have,  but  not 
a  necessity.  Tha  structure  outlined  in  this  research  effort  is  one 
possible  alternative  that  nay  warrant  additional  analysis.  Its  nain 
concern  has  bean  tha  patiant,  naasurad  by  tha  TIS  and  average  and 
naxinun  waiting  tines.  Thasa  perfornance  naasuras  ware  translated  into 
natarial  resources  required  to  achieve  those  desired  standards. 

Tha  need  for  a  casualty  evacuation  and  distribution  plan  cannot  be 
questioned,  even  one  that  is  sub-optinal.  A  raalistic,  workable  course 
of  action  affects  the  norale  of  tha  soldiers  and  their  support  fron  the 
people  at  hone.  The  will  to  fight  is  one  neasure  of  a  nation's  war¬ 
fighting  capability;  thus  actions  that  affect  that  willpower  clearly 
affect  the  outcone  of  tha  conflict. 
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Appendix  A 


SLAM  Network  Docunentation 

GEN,  WEUING,  THESIS,  10/22/85, ,  NO,  NO,  YES,  NO, ,  132) 
LIMITS, 10,2,2000> 

TIMST, XX( 4 ), CRAF  IN  USE, 8/10/1; 

TIMST,  XX(  3 ),  C?  IN  USE,10/0/l| 

TIMST,  NNQ(  2 ),  WAITING. 4. CRAF,  20/25/25; 

INTLC,  XX(  2  )«0 . 0,  XX<  4  >«0.0,  XX(  5  )-0.0,  XX(  4  )«12. 0; 


GLOBAL  VARIABLES 

XX(1)  -  ONE-DAY  CASUALTY  RATE(  X  .1) 

XX<  2  )  -  DAY  OF  WAR<  DOW ) 

XX<  3 )  ■  CASUALTY  COUNTER*  XX(  3 )  <  XX< 1 )  > 

XX(  4 )  -  NUMBER  OF  CRAF  IN  USE 
XX(  5 )  -  NUMBER  OF  C-9  IN  USE 
XX(  6 )  -  CRAF  CAPACITY* X  .1) 

XX<  7 )  -  STAGING  BASES!  3  ■  BWI,  4  -  NYC,  5  -  BOS,  4  »  PHL 
ATRIBUTES 

ATRIB(l)  -  START  TIME  FOR  PATIENT'S  TIME  IN  SYSTEM 
ATRIB<  2  >  -  PATIENT  CATEGORY 


NETWORK; 

RESOURCE/CRAF(  16 ),  8;  NUMBER  OF  B747  AIRCRAFT 

RES0URCE/C9(  9 ),  9;  NUMBER  OF  C-9  AIRCRAFT 

CREATE, 24,0,, 40,1;  40  ONE-DAY  TIME  PERIODS 

ASSIGN, XX<1)  -  TRIAG<  100.0, 175.0,250.0,1),!; 


VARIATIONS  FOR  GENERATING  CASUALTIES 

ASSIGN, XX(  1 )  *  UNFRM< 100.0,250.0, 1>,1; 
ASSIGN, XX( 1 )  -  RNORM* 173.0,25.0, 1  >,1; 
ASSIGN,  XX(  1 )  «  USERF(1),1; 


ASSIGN, XX(2)  *  XX(  2 )  +  1,XX(3)  ■  0.0,1; 

EV7  EVENT, 7,1;  DAILY  MTBED  UPDATE 

RTN  ASSIGN, XX< 3)  -  XX< 3  )  +  1,2; 

ACT , , XX( 3 )  .LT.  XX( 1  ),RTN;  GENERATES  DAILY  CASUALTIES 
ACT; 

QUEUE( 7  );  TEMPORARY  HOLDING  FILE 

ACT, 23.9/XX(  1 );  RELEASES  CASUALTIES  FOR  FLIGHT 

GOON, 1; 

ACT/2,,. 21, ASN1; CAT  A  CAS;  CAT  1  i  7  ARE  21%  OF  TOTAL 

ACT/3,,. 41, ASN2; CAT  B  CAS;  CAT  2,3,8  ARE  41%  OF  TOTAL 

ACT/4,,. 15, ASN3; CAT  C  CAS;  CAT  5,4,9,11  ARE  15%  OF  TOTAL 

ACT/5,, .23, ASN4; CAT  D  CAS;  CAT  4  1  10  ARE  23%  OF  TOTAL 

j - - - 

I 

;  THIS  SECTION  ASSIGNS  A  CATEGORY  TO  EACH  PATIENT 


ASN1  ASSIGN, ATRIBG 2)  -  1.0; 

ACT, ,,NXT; 

ASN2  ASSIGN, ATRIB<  2)  -  2.0; 

ACT, ,,NXT; 

ASN3  ASSIGN, ATRIB< 2)  -  3.0; 

ACT,,,NXT; 

ASN4  ASSIGN,  ATRIB(  2)  >  4.0; 

P . . 

;  THIS  SECTION  RELEASES  ONE  PLANE  LOAD  OF  PATIENTS  FOR  FLIGHT 
;  AND  CALLS  EVENT  1  TO  START  THE  CLOCK  FOR  TIME  IN  THE  SYSTEM. 

P 

NXT  GOON, 1; 

ACT , ,  NNGK  1 )  .CE.  XX(6),EV1; 

ACT,,, HOP; 

EV1  EVENT, 1,1; 

P . . 

;  THIS  SECTION  VERIFIES  THAT  A  LOAD  OF  PATIENTS  IS  READY 
;  AND  AT  LEAST  ONE  CRAF  AIRCRAFT  IS  AVAILABLE. 

P 

HOP  COON, 1; 

ACT , , NNQ( 2 >  .GE.  XX<6)  .AND.  NNRSC( 1 )  .CE.  1,QUE; 
ACT,,,Q1; 

QUE  QUEUE< 10 ), , , , SEL; 

P . . 

;  THIS  SECTION  ROUTES  EVERY  OTHER  CRAF  AIRCRAFT  TO  BUI 
P 

SEL  SELECT,, CYC,, QUE; 

ACT/20,,, PRI; ODD; 

ACT/1,,, SEC; EVEN; 

PRI  ASSIGN, XX( 7)  ■  3; 

ACT, , ,EV2; 

SEC  GOON, 1; 

ACT, , ,5,RGN4; 

ACT , , . 2, RGN5; 

ACT , , .3, RCN6; 

RGN4  ASSIGN,  XX(  7)  -  4; 

ACT, , ,EV2; 

RGN5  ASSIGN,  XX<  7)  -  5; 

ACT, , ,EV2; 

RGN4  ASSIGN,  XX<  7)  -  4; 

EV2  EVENT, 2,1; 

Q1  QUEUE<  1 ); 

TERM; 

P  - - 

;  FLY  CRAF  TO  FIRST  DESTINATION  HOSPITAL 


ENT1 


ENTER, 1,1; 

AUAIT( 8  ), CRAF/1, 1; 

ASSIGN, XX(4)  =  XX( 4 )  +  1.0; 

ACT, 3.5;  3.5 

GOON, 1; 

ACT/6, UNFRM< 8.5, 9.5,  2 >,ATRIB( 1  ) 
ACT/7, UNFRM(8. 0,9.0, 3), ATRIB<  1  > 
ACT/8, UNFRM(  7.7, 8.7, 4 ),ATRIB< 1 ) 
ACT/9, UNFRM< 8.2, 9.2, 5),  ATRIBG 1  ) 


HOUR  GROUND  TIME  BEFORE  T.O. 

.EQ.  3.0,BUI;FLY  TO  BUI; 

.EG).  4.0,NYC;FLY  TO  NYC; 

.EQ.  5.0,B0S;FLY  TO  BOS; 

.EQ.  6.0,PHL;FLY  TO  PHL; 


A  -  2 


BUI  EVENT, 3,1} 

TERMj 

NYC  EVENT, 4,1 1 

TERN} 

BOS  EVENT, 5,1 1 

TERMj 

PHL  EVENT, 6,1 j 

TERMj 

» . . . 

j  FLY  CRAF  TO  EUROPE  AND  FREE 

I 

ENT2  ENTER, 2,  lj 

ACVIO,  4.0  +  UNFRM(  7.0, 8.5,6  )jTO  EURj4  HR  GROUND  +  FLY  TIME 
FREE,CRAF/l,lj 
ASSIGN, XX(  4)  -  XX<4)  -  l.Oj 
TERMj 

I - - - - - 

j  FLY  CRAF  TO  I  FROM  DESTINATION  HOSPITAL 

j  THEN  RETURN  TO  EUROPE  AND  FREE 

I 

ENT3  ENTER, 3, lj 

ACT/11,1.5  +•  ATRIB(  1  )jCRAF  TO  H2j  1.5  HR  CC, REFUEL  +  FLY  TIME 
COON, lj 

ACT/12,3.0  +  ATRIB(1  )jCRAF  FM  H2j  3  HR  OFFLOAD  TIME  +  FLY  TIME 
GOON, 1; 

ACT/13, 4.0  +  UNFRM< 7.0, 8.5,7 >j TO  EURj4  HR  CC,REF,HX  ♦  FLY  TIME 
FREE, CRAF/1, 1 j 
ASSIGN, XX( 4)  »  XX(  4 )  -  l.Oj 
TERMj 

f - - - 

j  FLY  C9  TO  DESTINATION  HOSPITAL,  RETURN  AND  FREE 

I 

ENT4  ENTER, 4, lj 

AUAIT(  9  ),  C9/1, 1  j 

ASSIGN, XX< 5)  -  XX(  5 )  ♦  l.Oj 

ACT/14,2.3  +  ATRIB(  1  )j C9  TO  H2j  2.5  HR  ONLOAD  TIME  +  FLY  TIME 
GOON, lj 

ACT/15,2.0  ♦  ATRIB<  1  )jC9  FM  H2j  2  HR  OFFLOAD  TIME  +  FLY  TIME 
FREE,C9/1, lj 

ASSIGN, XX( 5)  -  XX< 5  )  -  l.Oj 
TERMj 

I  *  . . . . . 

j  COLLECT  STATS  FOR  PATIENTS  OFFLOADING  HOSP.l  VIA  CRAF 

I 

ENTS  ENTER, 3, lj 

COLCT,  INK  1 ),  TIS.H1  .CRAF, ,  1  j 

ACT,,,H1j 

I - 

j  COLLECT  STATS  FOR  OVERFLOW  PATIENTS 

J 

EN16  ENTER, 16,1; 

ACT/16,,, OVERjBWI  OVERFLOU; 

EN17  ENTER, 17,1; 

ACT/17,,, OVER;NYC  OVERFLOW; 


A  -  3 


EN18  ENTER, 18,  if 

ACT/18,,, OVERf BOS  OVERFLOW! 

EN19  ENTER, 19, If 

ACT/19,,, OVERf PHL  OVERFLOWf 
OVER  COLCT, INT< 1 ), OVERFLOW,,  If 

I  . . 

f  COLLECT  STATS  FOR  ALL  HOSP.l  PATIENTS 

I 

HI  COLCT, INK  1),TIS. HOSP.l,,  If 
ACT,,,AVCf 

I . . 

f  COLLECT  STATS  FOR  PATIENTS  OFFLOADING  H0SP.2+  VIA  CRAF/C9 

» 

ENT7  ENTER, 7, If 

ACT, 2.5  ♦  ATRIB(  2  )j 
COLCT, INK  1  ),TIS.H2.CRAF.C9, , If 
ACT, , ,H2f 

*  .  . . 

f  COLLECT  STATS  FOR  PATIENTS  OFFLOADING  H0SP.2+  VIA  CRAF 

I 

ENT8  ENTER, 8, If 

ACT, 1.5  *  ATRIB(  2  )f 
COLCT, INK  1 ),  TIS.H2.CRAF, , if 

f - - 

f  COLLECT  STATS  FOR  ALL  H0SP.2+  PATIENTS 

f 

H2  COLCT, INK  1 ), TIS.HOSP.2, ,  1; 

* . 

f  COLLECT  STATS  FOR  ALL  PATIENTS 

» 

AVC  COLCT,  INK  1),  AVERAGE,,  If 
TERMf 

ENDNETUORK; 

INIT, 0, 1 A40; 

SIMULATE; 

FINf 


Fortran  Docunantation 
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*  THESIS. FOR  -  10/22/85 
PROGRAM  MAIN 
DIMENSION  NSET< 20000) 

C0MM0N/SC0M1/ATRIB< 100 )f  DD(  100 ),  DDL(  100 ), DTNOU, II, UFA, MSTOP, NCLNR 
1, NCRDR, NPRNT, NNRUN, NNSET,  NTAPE,  SS(  100 >,SSL< 100),TNEXT,TN0U,XX< 100) 
COMMON  QSET( 20000) 

COMMON/MYVAR/MTBED< 3i 6, 30 >,  OUT(  3*6, 30, 61 ), FLYTMt 3»6, 30  ),FLY, IDEST 
EQUIVALENCE  ( NSETt 1 ), QSETt 1 ) ) 

NNSET  -  20000 

NCRDR-5 

NPRNT-6 

NTAPE-7 

NPLOT-2 

CALL  SLAM 

STOP 

END 

*  *  *  *  *  ttttttttttttttttttttttttt 

*  DEFINITION  OF  VARIABLES 

t 

%  MTBED(I,J)  ■  NUMBER  OF  EMPTY  BEDS  IN  HOSPITAL  J 

*  SERVED  BY  STAGING  BASE  I. 

*  OUT( I, J,K )  -  NUMBER  OF  PATIENTS  DISCHARGED  ON  DAY  K 

*  FROM  HOSPITAL  J  IN  STAGING  BASE  AREA  I. 

*  FLYTM(  I, J)  *  FLYING  TIME  FROM/TO  HOSPITAL  J 

*  TO/FROM  STAGING  BASE  I. 

*  FLY  *  DUMMY  VARIABLE  FOR  SELECTED  FLYING  TIME. 

*  IDEST  ■  DUMMY  VARIABLE  FOR  SELECTED  DESTINATION  HOSPITAL. 

*  AREA  •  ARRAY  SUBSCRIPT  FOR  STAGING  BASES. 

*  HOSP  -  ARRAY  SUBSCRIPT  FOR  HOSPITALS. 

*  DAY  -  ARRAY  SUBSCRIPT  FOR  DAY  OF  WAR. 

****************************** 

*  SUBROUTINE  TO  INITIALIZE  ALL  VARIABLES,  CONSTANTS  AND  ARRAYS. 

SUBROUTINE  INTLC 
DIMENSION  NSET< 20000) 

C0MM0N/SC0M1/ATRIB< 100 ),  DD< 100 ),  DDL< 100 ), DTNOU, 1 1, MFA, MSTOP, NCLNR 
1 , NCRDR, NPRNT, NNRUN,  NNSET,  NTAPE,  SS< 100  ), SSL< 100 ), TNEXT, TNOW, XXt 100  ) 
COMMON  QSET( 20000) 

EQUIVALENCE  ( NSET< 1 ),  QSET(  1  ) ) 

COMMON/MYVAR/MTBEDt  3 ! 6, 30  >,  OUT<  3 : 6, 30,.  61 ),  FLYTMt  3 1 6, 30  ),  FLY,  IDEST 
INTEGER  MTBED, OUT, AREA, HOSP, DAY 
REAL  FLY fM, FLY 
DO  120  AREA  -  3,6 
DO  110  HOSP  -  1,30 

MTBED( AREA,  HOSP)  =  0 
FLYTMt AREA, HOSP ?  *  0.0 
DO  100  DAY  -  1,61 

OUTt AREA, HOSP, DAY)  *  0 
100  CONTINUE 

110  CONTINUE 

120  CONTINUE 

* 


A  -  5 


.  *  * . •.'■j,'*)' .  •  ’«*'  .*•  .*«  ,** j,'-- »*• .  * 

Vi.*  frl/L-k'k.  V, 
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BALTIMORE,  WASHINGTON,  RICHMOND 
MTBED< 3, 1 )  •  436 

NORFQLK-VA  BEACH,  NEWPORT  NEWS-HAMPTON 
MTBED(3,2)  -  94 
FLYTM<  3, 2 )  ■  O.S 
PITTSBURGH 

MTBEDt 3, 3)  -  224 
FLYTMt  3,3)  -  0.6 
BUFFALO 

HTBED< 3,4)  -  118 
FLYTM<  3, 4 )  ■  0.8 
SYRACUSE 

HTBED( 3, 5  >  ■  42 
FLYTM(3,5)  -  0.8 
ROCHESTER 

HTBED<  3,6 )  -  56 
FLYTM< 3,6)  -  0.8 
ALBANY 

MTBED<  3,7 )  ■  74 
FLYTM<  3,7 )  ■  0.8 
RALEIGH-DURHAM 
MTBED( 3,8  )  -  61 
FLYTM( 3,8 )  -  0.9 
CLEVELAND 

MTBED( 3, 9 )  -  190 
FLYTMt 3, 9)  -  0.9 
COLUMBUS 

MTBEDt  3, 10 )  *  83 
FLYTMt 3, 10)  -  1.0 
DAYTON 

MTBEDt 3, 11)  ■  90 
FLYTMt 3, 11)  ■  1.1 
DETROIT, FLINT, TOLEDP 
MTBEDt 3, 12)  *  426 
FLYTMt 3, 12)  ■  1.2 
CINCINNATI 

MTBEDt 3, 13)  -  109 
FLYTMt 3, 13)  ■  1.2 
LEXINGTON-FAYETTE 
MTBEDt 3, 14)  -  44 
FLYTMt 3, 14)  *  1.2 
COLUMBIA 

MTBEDt 3, 15)  -  44 
FLYTMt 3, 15)  =  1.3 
KNOXVILLE 

MTBEDt 3, 16)  »  51 
FLYTMt 3, 16)  ■  1.3 
LOUISVILLE 

MTBEDt 3, 17)  «  88 
FLYTMt 3, 17)  *  1.4 
AUGUSTA 

MTBEDt 3, 18)  -  53 
FLYTMt 3, 18)  ■  1.4 


A  -  6 


*  INDIANAPOLIS 

MTBED( 3, 19 )  -  102 
FLYTM<3,19)  -  1.4 

*  GRAND  RAPIDS 

MTB£D<3,20)  -  42 
FLYTM<3,20)  -  1.5 

*  ATLANTA 

MTBED( 3,21 )  -  166 
FLYTM<3,21>  -  1.5 

*  NASHVILLE-DAVIDSON 

HTBED<  3, 22 )  -  99 
FLYTM<3,22>  -1.5 

*  CHICAGO, GARY 

MTBED<3,23>  ■  671 
FLYTM<  3,23 )  -  1.6 

*  MILWAUKEE 

MTBED( 3, 24  )  -  140 
FLYTM<3,24)  -  1.6 

*  BIRMINGHAM 

MTBED( 3, 25 )  -  93 
FLYTM<  3, 25 )  -  1.7 

*  JACKSONVILLE 

MTBED( 3, 26 )  -  55 
FLYTM( 3, 26  )  -  1.7 

*  ST  LOUTS 

MTbED( 3, 27  )  -  249 
FLYTM<  3, 27  )  -  1.8 

*  MEMPHIS 

MTBED( 3, 28 )  -  110 
FLYTM<  3, 28 )  -  1.9 

*  TAMPA, ORLANDO 

MTBED( 3, 29 )  -  206 
FLYTM(  3, 29 )  *  2.3 
♦TOTAL  BUI  AREA  »  4216 

* 

*  NEW  YORK,  NEW  HAVEN-UEST  HAVEN,  JERSEY  CITY,  NEU  BRUNSWICK, 

*  LONG  BRANCH-ASHBURY  PARK,  PATTERSON-CLIFTON,  NEWARK, 

*  POUGHKEEPSIE,  NEWBURGH-MIDDLETOWN,  NASSAU-SUFFOLK 

MTBED<  4, 1 >  ■  1342 

*  MINNEAPOLIS-ST  PAUL 

MTBED( 4,3 )  -  186 
FLYTM( 4,3 )  -  2.0 

*  DES  MOINES 

MTBED<  4, 4  )  •  40 
FLYTM( 4, 4 )  »  2.1 
%  OMAHA 

MTBED( 4,5)  *  75 
FLYTH(  4,5)  ■  2.3 

*  TOTAL  NYC  AREA  =  1643 
$ 

%  BOSTON,  BROCKTON,  HARTFORD,  NEW  LONDON-NORWICH,  PROVIDENCE 

MTBED( 5, 1 )  -  441 

*  KANSAS  CITY,  TOPEKA 

MTBED<  5, 4  )  *  137 
FLYTM(  5, 4  )  >  2.4 
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*  TULSA 

MTBEDt  5, 5)  -  55 
FLYTMt 5, 5 )  ■  2.7 

*  WICHITA 

MTBEDt  5, 6)  -  44 
FLYTM(  5, 6 )  ■  2.8 

*  OKLAHOMA  CITY 

MTBED( 5,7 )  ■  129 
FLYTH<  5, 7 )  -  3.0 

*  TOTAL  BOS  AREA  -  826 

*  PHILADELPHIA,  WILMINGTON,  TRENTON 

MTBED<  6, 1 )  -  439 

*  MOBILE,  BILOXI -GULFPORT 

MTBEDt  6, 3)  *  51 
FLYTM( 6, 3 >  •  1.9 
t  JACKSON 

MTBEDt  6, 4)  -  45 
FLYTM<  6, 4 )  -  1.9 

*  LITTLE  ROCK 

MTBED<  6, 5  >  -  68 
FLYTM<  6, 5 )  ■  2.0 

*  HI AMI,  FT  LAUDERDALE 

MTBEDt 6, 6)  -  257 
FLYTMt  6,6 )  -  2.0 

*  NEW  ORLEANS 

MTBED<6,7)  -  123 
FLYTMt 6, 7)  -  2.2 

*  SHREVEPORT 

HTBEDt  6,8 )  -  49 
FLYTMt 6, 8)  -  2.3 

*  DALLAS 

MTBEDt  6, 9 )  -  219 
FLYTMt 6, 9)  -2.5 

*  HOUSTON 

MTBEDt 6, 10)  «  244 
FLYTMt 6, 10)  ■  2.6 

*  AUSTIN 

MTBEDt 6, 11)  ■  26 
FLYTMt 6, 11)  ■  2.8 

*  SAN  ANTONIO 

MTBEDt 6, 12)  »  112 
FLYTMt 6, 12)  «  2.9 

*  TOTAL  PHL  AREA  =  1633 

*  GRAND  TOTAL  ■  8318 

RETURN 

END 

* 
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*************** 

*  MORE  VARIABLES  DEFINED 

* 

*  DOW  -  DAY  OF  WAR  (EQUIVALENCE  XX<2>). 

*  CITY  -  STACINC  BASE  (EQUIVALENCE  XX(7)). 

*  CAPY  -  CAPACITY  OF  CRAFt X  .1) 

*  RCN  -  REGION  FOR  EACH  STAGING  BASE. 

*  NRSC  -  NUMBER  OF  RESOURCES(  SLAM  VAR). 

*  TMP  •  DUMMY  VARIABLE. 

*  CAP  ■  C-9  CAPACITY(X  .1) 

****************************** 

*  EVENT  SUBROUTINES 

* 

SUBROUTINE  EVENT(N) 

DIMENSION  NSET(  20000) 

COMMON/SCOHl/ATRIBt  100  >,  DD(  100  >,  DDL(  100 ),  DTNOW,  I I, MFA,  MSTOP,  NCLNR 
1,  NCRDR, NPRNT, NNRUN,  NNSET, NTAPE, SS( 100 ), SSL( 100 ),  TNEXT,  TNOW,  XX(  100  ) 
COMMON  QSET( 20000) 

COMMON/MYVAR/MTBED(  3*6, 30 ),OUT( 3 i 6, 30,41  ),FLYTM(  3>6, 30 )tFLYf IDEST 
EQUIVALENCE  ( NSET(  1 ), QSET(  1 ) ) 

EQUIVALENCE  (XX(2),D0W> 

EQUIVALENCE  (XXU),CAPY) 

EQUIVALENCE  (XX(7),CITY> 

INTEGER  MTBED,  OUT,  AREA,  HOSP,  DAY,  CAPY,  RCN, NRSC,  TMP,  MTBD,  CAP 
REAL  FLYTM, FLY,  CITY 
CAP  -  4 

GOTO  ( 1, 2,3, 4,5, 4,7 )  N 

*************** 

*  EVENT  1 

* 

*  WHEN  FULL  CRAF  LOAD  IS  AVAILABLE,  REMOVE  PATIENTS  FROM  FILE  1. 

* 

1  DO  130  IA  -  CAPY, 1,-1 

CALL  RMOVE( IA, 1, ATRIB  ) 

ATRIB(l)  ■  TNOW 
CALL  FILEM( 2, ATRIB ) 

130  CONTINUE 
RETURN 

*************** 

*  EVENT  2 

*  DETERMINE  DESTINATION  REGION, 

*  REMOVE  PATIENTS  FROM  FILE  2  FOR  FLIGHT, 

*  FLY  CRAF  TO  U.S. 

*  FILE  3  *  BUI, SOX  I  FILE  4  -  NYC,25X 

*  FILE  3  -  BOS, 10X  I  FILE  6  -  PHL,15X 

2  RGN  -  INT( CITY  ) 

DO  140  IB  »  CAPY, 1,-1 
CALL  RMQVE(  IB, 2,  ATRIB ) 

CALL  FILEM( RGN, ATRIB) 

140  CONTINUE 


**********  *  *  *  * 

EVENT  3 


*  PATIENTS  AND  CRAF  HAVE  ARRIVED  AT  BUI 

*  FLY  CRAF  TO  EUROPE  AND  FREE. 

*  IF  AT  LEAST  752  OF  ONE  C-9  LOAD  ARE  IN  FILE  3  AND 

*  IF  "SERCH"  RETURNS  WITH  A  DESTINATION  HOSPITAL  AND 

*  IF  AT  LEAST  1  C-9  IS  AVAILABLE  THEN 

*  REMOVE  PATIENTS  FROM  FILE  3  (UPLOAD  AIRCRAFT) 

*  UPDATE  MTBED  AND  DETERMINE  DAY  OF  DISCHARGE, 

*  FLY  PATIENTS  TO  DESTINATION  AND  COLLECT  STATS, 

*  SCHEDULE  C-9  FLIGHT, 

*  CHECK  FOR  MORE  PATIENTS,  ETC) 

*  ELSE  OFFLOAD  ALL  REMAINING  PATIENTS, 

*  ELSE  OFFLOAD  ALL  REMAINING  PATIENTS, 

*  ELSE  OFFLOAD  ALL  REMAINING  PATIENTS. 

3  CALL  ENTER( 2, ATRIB ) 

NRSC  -  NNRSCC2) 

TMP  -  CAPY 

150  IF  (TMP  .GE.  3)  THEN 

IF  (TMP  .Ed.  3)  THEN 
CAP  ■  3 
ENDIF 

CALL  SERCH(  3, CAP ) 

IF  (IDEST  .GE.  2)  THEN 
IF  (NRSC  .GE.  1)  THEN 
HOSP  -  IDEST 
DO  160  IC  «  CAP, 1,-1 

CALL  RMQVE(  IC, 3,  ATRIB ) 

CALL  CASE( 3, HOSP  ) 

ATRIB( 2 )  ■  FLY 
CALL  ENTER(  7,  ATRIB  > 

160  CONTINUE 

NRSC  -  NRSC  -  1 
TMP  ■  TMP  -  CAP 
ATRIB(l)  -  FLY 
CALL  ENTER( 4, ATRIB  ) 

GOTO  150 
ELSE 

CALL  OFLD<  3, TMP ) 

ENDIF 
ELSE 

CALL  QFLD( 3, TMP ) 

ENDIF 
ELSE 

CALL  OFLD(  3, TMP ) 

ENDIF 
RETURN 

*************** 

*  EVENT  4,  5,  OR  6 

* 

*  PATIENTS  AND  CRAF  HAVE  ARRIVED  AT  NYC,  BOS  OR  PHL. 

* 
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*  IF  AT  LEAST  "CAPY"  EMPTY  BEDS  ARE  AVAILABLE  THEN 

*  RETURN  CRAF,  OFFLOAD  PATIENTS  AND  GET  STATS, 

*  ELSE  FIND  FIRST  HOSPITAL  UITH  ENOUGH  MTBEDS, 

*  IF  ONE  IS  FOUND  THEN 

*  SCHEDULE  CRAF  TO  l  FROM  BEST  HOSP,  RETURN  TO  EUROPE, 

*  REMOVE  PATIENTS  FROM  FILE  "AREA", 

*  UPDATE  MTBEDS  AND  DAY  OF  DISCHARGE, 

*  FLY  PATIENTS  TO  HOSPITAL,  GET  STATS. 

*  ELSE  OFFLOAD  PATIENTS  (POTENTIAL  OVERFLOW) 

* 

4  AREA  -  4 

GOTO  170 

5  AREA  -  5 

GOTO  170 

6  AREA  -  6 

* 

170  MTBD  -  MTBED(  AREA,  1 ) 

IF  (MTBD  .GE.  CAPY)  THEN 
CALL  ENTER( 2, ATR1B ) 

CALL  OFLD( AREA, CAPY) 

ELSE 

CALL  SERCH( AREA, CAPY) 

IF  (IDEST  .GE.  2)  THEN 
HOSP  -  IDEST 
ATRIB(l)  =  FLY 
CALL  ENTER( 3, ATRIB  ) 

DO  180  ID  -  CAPY, 1,-1 

CALL  RMOVE(  ID, AREA, ATRIB) 

CALL  CASE(  AREA,  HOSP) 

ATRIB(2)  -  FLY 
CALL  ENTER(  8, ATRIB  ) 

180  CONTINUE 

ELSE 

CALL  OFLD( AREA, CAPY) 

ENDIF 

ENDIF 

RETURN 

*************** 

*  EVENT  7 

* 

*  DAILY  MTBED  UPDATE 

* 

*  MTBED  =  EMPTY  BEDS  AT  HOSPITAL  DEFINED  BY  (AREA, HOSP) 

*  OUT  *  NUMBER  OF  PATIENTS  DISCHARGED  FROM  HOSPITAL  ON  A 

*  SPECIFIED  DAY  OF  THE  WAR  (DOW) 

* 

7  DAY  *  INT(DOW) 

DO  200  AREA  >3,4 

DO  190  HOSP  =1,30 

MTBED(  AREA,  HOSP  )  =  MTBED(  AREA,  HOSP  )  +  OUT(  AREA, HOSP,  DAY  ) 
190  CONTINUE 

200  CONTINUE 
RETURN 
END 

****************************** 
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*  SUBROUTINE  OFLD  REMOVES  "NOFF"  PATIENTS  FROM  FILE  "AREA",  UPDATES 

*  MTBEDS  AND  COLLECTS  STATS  FOR  ALL  PATIENTS  INCLUDING  OVERFLOW. 

* 

SUBROUTINE  OFLD(  AREA, NOFF  > 

DIMENSION  NSET( 20000) 

C0MM0N/SCQM1/ATRIB<  100  >,  DD<  100  >,  DDL(  100 ),  DTNOU,  II,  MFA,  HSTOP,  NCLNR 
1 ,  NCRDR,  NPRNT,  NNRUN,  NNSET,  NTAPE,  SS(  100  ),  SSL(  100 ),  TNEXT,  TNOU,  XX<  100  ) 
COMMON  QSET( 20000) 

COMMON/MYVAR/MTBED(  3«6,30 ), OUT<  3 1 6, 30, 61 ), FLYTM(  3 1 6, 30 ), FLY,  IDEST 
EQUIVALENCE  ( NSET< 1 ),QSET(  1 ) ) 

EQUIVALENCE  (XX(2),D0U) 

INTEGER  MTBED,  OUT,  AREA,  HOSP,  DAY,  NOFF,  OVER,  MTB 
IF  (NOFF  .LE.  0)  THEN 
RETURN 
ENDIF 

MTB  -  MTBED< AREA,  1 ) 

* 

* 

IF  (MTB  .GE.  NOFF)  THEN 
DO  210  JA  -  NOFF, 1,-1 

CALL  RMOVE(  JA, AREA,  ATRIB  > 

CALL  CASE( AREA, 1 ) 

CALL  ENTER( 5, ATRIB ) 

210  CONTINUE 

ELSE 

OVER  -  NOFF  -  MTB 
DO  220  JB  -  MTB, 1,-1 

CALL  RMOVE(  JB,  AREA,  ATRIB  ) 

CALL  CASE( AREA, 1 ) 

CALL  ENTER(3,ATRIB) 

220  CONTINUE 

DO  230  JC  -  OVER, 1,-1 

CALL  RMOVE(  JC,  AREA,  ATRIB ) 

CALL  ENTER( AREA  ♦  13, ATRIB) 

230  CONTINUE 

ENDIF 
RETURN 
END 

*****«*t*t****t**»tt******t*tt 

*  SUBROUTINE  CASE  DECREASES  THE  NUMBER  OF  MTBEDS  AND  INCREASES  THE 

»  NUMBER  OF  PATIENTS  TO  BE  DISCHARGED  ON  A  PARTICULAR  DAY  OF 

*  THE  UAR  (DOW)  DEPENDING  ON  THEIR  MEDICAL  CATEGORY  ( ATRIB( 2 ) ) 

SUBROUTINE  CASEt  AREA, HOSP ) 

DIMENSION  NSET< 20000) 

C0MH0N/SC0M1/ATRIB(  100 ),  DD(  100  ),  DDL(  100  ),  DTNOW,  II, MFA,  MSTOP,  NCLNR 
1,  NCRDR,  NPRNT,  NNRUN,  NNSET,  NTAPE,  SS<  100  ),  SSL(  100  ),  TNEXT,  TNOU,  XX(  100  ) 
COMMON  QSET(  20000) 

COHMON/MYVAR/MTBED(  3 » 6, 30  ),  OUT(  3 » 6, 30, 61  >,  FLYTM(  3:6,30),  FLY,  IDEST 
EQUIVALENCE  ( NSET< 1 ), QSET< 1 ) ) 

EQUIVALENCE  <XX(2),D0W) 

INTEGER  MTBED, OUT, AREA, HOSP, DAY 

t 

i  UPDATE  MTBED  ARRAY 

MTBED( AREA, HOSP)  -  MTBED(  AREA, HOSP  )  -  1 
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IF  ( ATRIBf 2 )  .EQ.  1.0)  THEN 
DAY  INK  DOW  ♦  16) 

IF  <  DAY  .LE.  60)  THEN 

OUT(  AREA,  HOSP,  DAY )  -  OUT<  AREA,  HOSP,  DAY )  +  1 
ELSE 

OUT<  AREA,  HOSP,  61)  -  OUT(  AREA,  HOSP,  61 )  ♦  1 
ENDIF 

ELSEIF  ( ATRIB< 2 )  .EQ.  2.0)  THEN 
DAY  ■  INK  DON  ♦  25) 

IF  (DAY  .LE.  60)  THEN 

OUT( AREA, HOSP, DAY)  -  OUT ( AREA, HOSP, DAY )  ♦  1 
ELSE 

OUT<  AREA,  HOSP,  61)  -  OUT<  AREA,  HOSP,  61 )  ♦  1 
ENDIF 

ELSEIF  ( ATRIB< 2 )  .EQ.  3.0)  THEN 
DAY  •  INK  DOW  ♦  38) 

IF  (DAY  .CT.  60)  THEN 

OUT(  AREA,  HOSP,  61)  -  OUT(  AREA,  HOSP,  61 )  +  1 
ELSE 

OUT(  AREA,  HOSP,  DAY)  -  OUT(  AREA,  HOSP,  DAY )  +  1 
ENDIF 
ELSE 

DAY  -  INK  DOW  +  51) 

IF  (DAY  .CT.  60)  THEN 

OUT(  AREA,  HOSP,  61)  -  OUT(  AREA,  HOSP,  61 )  +  1 
ELSE 

OUT(  AREA,  HOSP,  DAY)  -  OUT(  AREA,  HOSP,  DAY  )  ♦  1 
ENDIF 
ENDIF 
RETURN 
END 

*t**tt**t*t*tt***t*tt*t***t** 

USER  FUNCTION  ALLOWS  CASUALTIES  TO  BE  GENERATED  IN  A  CYCLIC 
PATTERN.  THIS  EXAMPLE  SHOW  A  10  DAY  CYCLE  WHERE  THE 
QUANTITY  2*PI*10*D0W  -  ARGUMENT,  175  -  EXPECTED  VALUE  AND 
175  CORRESPONDS  TO  THE  MAX  DEVIATION  FROM  EXPECTED  VALUE. 


FUNCTION  USERF(I) 

DIMENSION  NSET(  20000) 

COMMON/SCOM1/ATRIB<  100),DD<  100>,DDL(  100  ),  DTNOU,  II,MFA,MSTOP,NCLNR 
1,  NCRDR,  NPRNT,  NNRUN,  NNSET,  NTAPE,  SS(  100  ),SSL(  100),TNEXT,TNOU,XX(  100  ) 
COMMON  QSET(  20000) 

COMMON/MYVAR/MTBEDC  3 1 6, 30  ), OUT(  3! 6,  30,  61  ),  FLYTM(  3 S6, 30 ), FLY,  IDEST 
EQUIVALENCE  ( NSET< 1 ), OSETC 1 > ) 

EQUIVALENCE  <XX(2),D0W> 

GOTO  (1)1 

1  USERF  »  175  ♦  175*SIN(0.6283*DOW>  / 

RETURN  / 

END  ' 


t 

% 

t 

% 

t 


ft**************************** 

SUBROUTINE  SERCH  FINDS  THE  NEAREST  HOSPITAL  WITH  AT  LEAST  "BEDS 
NUMBER  OF  BEDS  AVAILABLE  WITHIN  EACH  AREA.  RETURNS  WITH 


HOSPITAL  ID  "IDEST"  AND  FLYTIME  "FLY". 
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SUBROUTINE  SERCH< AREA, BEDS ) 

DIMENSION  NSETt 20000) 

COMMON/SCOM1/ATRIB<  100)tDD<  100),  DDK  100  ),  DTNOU,  1 1 ,  MFA,  MSTOP,  NCLNR 
1,  NCRDR,  NPRNT,  NNRUN,  NNSET,  NTAPE,  SS(  100 ),  SSL<  100  ),  TNEXT,  TNOH,  XX(  100 ) 
COMMON  QSET( 20000) 

COMMON/M YVAR/HTBEIX  316,30),  OUTi  3 * 6 , 30 , 6 1  )f  FLYTM(  3 « 6, 30), FLY,  IDEST 
EQUIVALENCE  ( NSETt 1 ),QS£T< 1 )  > 

INTEGER  MTBED, AREA, HOSP,  BEDS 
REAL  FLYTM,FLY 

* 

IDEST  -  1 
FLY  ■  0.0 
DO  240  HOSP  >  2,30 

IF  <HTBED< AREA, HOSP)  .GE.  BEDS)  THEN 
IDEST  «  HOSP 
FLY  -  FLYTM( AREA, HOSP) 

RETURN 
ENDIF 
240  CONTINUE 
RETURN 
END 

****************************** 

SUBROUTINE  OTPUT 
DIMENSION  NSET( 20000) 

COMMON/SCOMl/ATRIBt  100 ),  DD<  100  ),  DDL(  100  >,  DTNOU,  II ,  MFA,  MSTOP,  NCLNR 
1,  NCRDR,  NPRNT,  NNRUN,  NNSET,  NTAPE,  SS<  100  ),  SSL<  100 ),  TNEXT,  TNOU,  XX(  100  ) 
COMMON  QSET( 20000) 

COMMON/MYVAR/MTBED(  3 30 ),  OUT<  3*6,30,61),  FLYTh(  3 *6, 30 ),  FLY,  IDEST 
EQUIVALENCE  ( NSET(  1 ), QSET(  1 ) ) 

INTEGER  MTBED, OUT 
REAL  FLYTM,FLY 
RETURN 
END 


APPENDIX  B 


POTENTIAL  NDHS  PARTICIPANTS 


cm 

81 

BEDS 

1. 

Mobil* 

AL 

1325 

2. 

Birninghan 

AL 

3734 

3. 

Littl*  Rock 

AR 

2732 

4. 

Phoenix 

AZ 

4100 

5. 

Tuscon 

AZ 

1659 

6. 

Riverside 

CA 

3758 

7. 

San  Jose 

CA 

3353 

8. 

Sacranento 

CA 

1708 

9. 

Stockton 

CA 

654 

10. 

Santa  Rosa 

CA 

542 

11. 

Los  Ang*l*s-Long  Beach 

CA 

22600 

12. 

San  Francisco-Qakland 

CA 

9273 

13. 

San  Diego 

CA 

4766 

14. 

Anahein-Santa  Ana 

CA 

4430 

13. 

Denver 

CO 

5283 

16. 

New  Haven-West  Haven 

CT 

2443 

17. 

Hartford 

CT 

2332 

18. 

New  London-Norwich 

CT 

421 

19. 

Washington 

DC 

8119 

20. 

UilMington 

DE 

1544 

21. 

Ft  Lauderdale 

FL 

3741 

22. 

Orlando 

FL 

2149 

23. 

Jacksonville 

FL 

2239 

24. 

Miani 

FL 

6359 

23. 

Tanpa-St  Petersburg 

FL 

6119 

26. 

Augusta 

GA 

2144 

27. 

Atlanta 

GA 

6657 

28. 

Des  Moines 

IA 

1626 

29. 

Chicago 

IL 

24866 

30. 

Gary-Hannond 

IN 

1979 

31. 

Indianapolis 

IN 

4104 

32. 

Wichita 

KS 

1769 

33. 

Topeka 

KS 

677 

34. 

Louisville 

KY 

3559 

33. 

Lexington-Fayette 

KY 

1784 

36. 

New  Orleans 

LA 

4924 

37. 

Shreveport 

LA 

1973 

38. 

Brockton 

MA 

576 

39. 

Boston 

MA 

12086 

40. 

Baltinore 

MD 

6447 

41. 

Flint 

MI 

1583 

42. 

Grand  Rapids 

MI 

1682 
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43. 

Detroit 

HZ 

12968 

44. 

Hinneapolis-St  Paul 

KN 

7440 

45. 

St  Louis 

HO 

9992 

46. 

Kansas  City 

HO 

5628 

47. 

Biloxi-Gulfport 

H8 

754 

48. 

Jackson 

H8 

1819 

49. 

Raleigh-Durhan 

NC 

2474 

50. 

Onaha 

NE 

3010 

51. 

Jarsay  City 

NJ 

1642 

52. 

Naw  Brunswick 

NJ 

1321 

53. 

Long  Branch-Asbury  Park 

NJ 

1208 

54. 

Tranton 

NJ 

1184 

55. 

Pattarson-Clifton 

NJ 

1177 

56. 

Nawark 

NJ 

6689 

57. 

Albuquerque 

NH 

1508 

58. 

Las  Vegas 

NV 

1174 

59. 

Poughkaapsia 

NY 

777 

60. 

Nawburgh-Hiddlaton 

NY 

772 

61. 

Buffalo 

NY 

4741 

62. 

Albany-Schenactady 

NY 

2973 

63. 

Rochester 

NY 

2258 

64. 

Syracusa 

NY 

1701 

65. 

Naw  York 

NY 

31331 

66. 

Nassau-Suffolk 

NY 

6353 

67. 

Cincinnati 

OH 

4363 

68. 

Dayton 

OH 

3635 

69. 

Colunbus 

OH 

3342 

70. 

Tolado 

OH 

2522 

71. 

Clavaland 

OH 

7622 

72. 

Oklahona  City 

OK 

3535 

73. 

Tulsa 

OK 

2223 

74. 

Portland 

OR 

3688 

75. 

Philadelphia 

PA 

14863 

76. 

Pittsburgh 

PA 

8969 

77. 

Providanca 

RI 

2252 

78. 

Coluhbi* 

SC 

1767 

79. 

Metiphis 

TN 

4411 

80. 

Nashville-Davidson 

TN 

3961 

81. 

Knoxvilla 

TN 

2076 

82. 

San  Antonio 

TX 

4505 

83. 

El  Paso 

TX 

1647 

84. 

Austin 

TX 

1062 

85. 

Houston 

TX 

9790 

86. 

Dallas-Ft  Worth 

TX 

8785 

87. 

Salt  Lake  City 

UT 

2205 

88. 

Richnond 

VA 

2878 

89. 

Nor folk -Va  Beach 

VA 

2516 

90. 

Newport  News-Hanpton 

VA 

1244 

91. 

Seattle-Everett 

UA 

2998 

92. 

Tacona 

UA 

1134 

93. 

Spokane 

UA 

1070 

94. 

Milwaukee 

WI 

5628 
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APPENDIX  C 


4  Factor  Analysis  of  Variance  for  Mean  Tine  in  Systen  (2  replicates). 


Source 

Sun  of 

Degrees  of 

Mean 

F 

Tail 

Squares 

Freedon 

Square 

Stat 

Prob 

Cas 

2.41501 

1 

2.41501 

179.24 

.00 

Craf 

1.25696 

1 

1.25696 

93.29 

.00 

C9 

0.02838 

1 

0.02838 

2.11 

.16 

Cap 

2.84214 

1 

2.84214 

210.94 

.00 

Cas/Craf 

0.53510 

1 

0.53510 

39.71 

.00 

Cas/C9 

0.00455 

1 

0.00455 

.34 

.57 

Craf/C9 

0.00617 

1 

0.00617 

.46 

.51 

Cas/Cap 

1.92624 

1 

1.92624 

142.97 

.00 

Craf /Cap 

0.84737 

1 

0.84737 

62.89 

.00 

C9/Cap 

0.00948 

1 

0.00948 

.70 

.41 

Cas/Craf/C9 

0.01464 

1 

0.01464 

1.09 

.31 

Cas/Craf /Cap 

0.30890 

1 

0.30890 

22.93 

.00 

Cas/C9/Cap 

0.00525 

1 

0.00525 

.39 

.54 

Craf/C9/Cap 

0.00316 

1 

0.00316 

.23 

.63 

Cas/Craf /C9/Cap  0.01867 

1 

0.01867 

1.39 

.26 

Error 

0.21557 

16 

0.21557 

3  Factor  Analysis  of  Variance  for  Mean  Tine  in  Systen  (3  replicates). 


Source 

Sun  of 
Squares 

Degrees 

Freedon 

of  Mean 

Square 

F 

Stat 

Tail 

Prob 

Cas 

2.38442 

1 

2.38442 

365.21 

.00 

Craf 

0,77676 

1 

0.77676 

118.97 

.00 

Cap 

2.39895 

1 

2.39895 

367.43 

.00 

Cas/Craf 

0.31116 

1 

0.31116 

47.66 

.00 

Cas/Cap 

1.73454 

1 

1.73454 

265.67 

.00 

Craf /Cap 

0.41582 

1 

0.41582 

63.69 

.00 

Cas/Craf/Cap 

0.11174 

1 

0.11174 

17.11 

.00 

Error 

0.10446 

16 

0.00653 

;  / 


Regression  Statistics  for  Three-Factor  Linear  Model 


Adjusted  Squared  Mult.  Corr. 
Residual  Mean  Square 

Standard  Error  of  Est. 
F-Statistic 

Nunerator  Degrees  of  Freedon 
Denoninator  Degrees  of  Freedon 
Significance  (Tail  prob. ) 

.98535 

.00519 

.07204 

221.96 

7 

16 

.00000 

Variable 

Regression 

Stand. 

r- 

Contribution 

Nane 

Coefficient 

Coeff. 

Stat 

to  R-Square 

Intercept 

3.03254 

5.095 

206.22 

Cas 

0.309242 

0.531 

21.02 

.282 

Craf 

-0.179915 

-0.309 

-12.23 

.095 

Cas/Craf 

-0.108445 

-0.186 

-  7.37 

.035 

Cap 

-0.325440 

-0.559 

-22.14 

.312 

Cas/Cap 

-0.241037 

-0.448 

-17.75 

.201 

Craf/Cap 

0.131528 

0.226 

8.94 

.051 

Cas/Craf /Cap 

0.070494 

.121 

4.49 

.015 

Regression  Statistics  for  Linear 

Model  (16 

aircraft ). 

Adjusted  Squared  Mult.  Corr. 

.98619 

Residual  Mean  Square 

.00758 

Standard  Error  of  Est. 

.08707 

F-Statistic 

262.84 

Nunerator  Degrees  of  Freedon 

3 

Denoninator  1 

Degrees  of  Freedon 

8 

Significance 

<  Tail  Prob. ) 

.00000 

Variable 

Regression 

Stand. 

T- 

Contribution 

Nane 

Coefficient 

Coeff . 

Stat 

to  R-Square 

Intercept 

3.21243 

4.335 

127.80 

Cas 

0.429043 

0,605 

17.07 

.366 

Cap 

-0.447786 

-0,631 

-17.81 

.398 

Cas/Cap 

-0.337070 

-0.475 

-13.41 

.226 

C  -  2 
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Regression  Statistics  for  Linear  Model  (17  aircraft) 


Adjusted  Squared  Mult.  Corr.  .99421 
Residual  Mean  Square  .00204 
Standard  Error  of  Est.  .04524 
F-Statistic  630.28 
Nunerator  Degrees  of  Freedon  3 
Denoninator  Degrees  of  Freedon  8 
Significance  ( Tail  Prob. )  .00000 


Variable 

Regression 

Stand. 

T- 

Contribution 

bane 

Coefficient 

Coeff. 

Stat 

to  R-Square 

Intercept 

3.02838 

5.095 

231.90 

Cas 

0.336020 

0.590 

25.73 

.349 

Cap 

-0.342767 

-0.602 

-26.25 

.363 

Cas/Cap 

-0.303418 

-0.533 

-23.23 

.284 

Regression 

Statistics  for  Linear 

Model  ( 18 

aircraft ). 

Adjusted  Squared  Mult.  Corr.  .96679 
Residual  Mean  Square  .00437 
Standard  Error  of  Est.  .06614 
F-Statistic  107.73 
Nunerator  Degrees  of  Freedon  3 
Denoninator  Degrees  of  Freedon  8 
Significance  (Tail  Prob.)  .00000 


Variable 

Regression 

Stand. 

T- 

Contribution 

Nane 

Coefficient 

Coeff. 

Stat 

to  R-Square 

Intercept 

2.  F, 6087 

7.883 

149.85 

Cas 

0.202296 

0.582 

10.60 

.339 

Cap 

-0.200557 

-0.577 

-10.50 

.333 

Cas/Cap 

-0.191466 

-0.551 

-10.03 

.304 

Regression  Statistics  for  Three-Factor  Quadratic  Hodel 


Adjusted  Squared  Mult,  Corr.  .95046 
Residual  Mean  Square  .01125 
Standard  Error  of  Est.  .10606 
F-Statistic  220.27 
Numerator  Degrees  of  Freedon  7 
Denoninator  Degrees  of  Freedon  73 
Significance  <Tail  prob. )  .00000 


Variable 

Regression 

Stand. 

T-  Contribution 

Name 

Coefficient 

Coeff. 

Stat  to  R-Square 

Intercept 

-2.979860 

5.521 

-  .54 

Cas 

2.095410 

5.419 

4.78 

.014 

Cap 

-1.810300 

-3.121 

-  2.71 

.005 

Cas/Cas 

0.061310 

5.551 

5.52 

.019 

Cap/Cap 

0.096583 

3.998 

3.86 

.009 

Cas/Craf 

-0.107817 

-5.754 

-12.75 

.101 

Cas/Cap 

-0.183286 

-7.943 

-15.55 

.150 

Craf/Cap 

0.140813 

5.056 

11.41 

.081 

Regression  Statistics  for  Quadratic  Model  (16  aircraft). 


Adjusted  Squared  Mult,  Corr.  .95983 
Residual  Mean  Square  .01476 
Standard  Error  of  Est.  .12151 
F-Statistic  208.08 
Nunerator  Degrees  of  FreedoM  3 
Denoninator  Degrees  of  Freedom  23 
Significance  (Tail  Prob.)  .00000 


Variable 

Regression 

Stand. 

T-  Contribution 

Name 

Coefficient 

Coeff. 

Stat  to  R-Square 

Intercept 

3.100320 

5.114 

12.78 

Cas/Cas 

0.081006 

5.839 

11.27 

.196 

Cap/Cap 

0.134249 

4.423 

8.78 

.119 

Cas/Cap 

-0.210844 

-7.273 

-10.08 

.157 

Regression  Statistics  for  Quadratic  Model  (17  aircraft) 


R 


Adjusted  Squared  Mult.  Corr. 
Residual  Mean  Square 

Standard  Error  of  Est. 
F-Statistic 

Numerator  Degrees  of  Freedom 
Denominator  Degrees  of  Freedom 
Significance  (Tail  Prob. ) 

.97187 

.00522 

.07226 

300.43 

3 

23 

.00000 

Variable 

Regression 

Stand. 

T- 

Contribution 

Name 

Coefficient 

Coeff. 

Stat 

to  R-Square 

Intercept 

Cas/Cas 

2.751400 

0.073736 

6.386 

7.478 

19.07 

17.25 

.322 

Cap/Cap 

0.131641 

6.103 

14.48 

.227 

Cas/Cap 

-0.197806 

-9.601 

-15.90 

.274 

Regression  Statistics  for  Quadratic  Model  (18  aircraft). 


Adjusted  Squared  Mult.  Corr.  .86693 
Residual  Mean  Square  .00845 
Standard  Error  of  Est.  .09190 
F-Statistic  57.46 
Numerator  Degrees  of  Freedom  3 
Denoninator  Decrees  of  Freedom  23 
Significance  (  fail  Prob. )  .00000 


Variable 

Regression 

Stand. 

T-  Contribution 

Name 

Coefficient 

Coeff. 

Stat  to  R-Square 

Intercept 

2.545690 

10.105 

13.88 

Cas/Cas 

0.046061 

7.989 

8.47 

.368 

Cap/Cap 

0.084815 

6.725 

7.34 

.275 

Cas/Cap 

-0.124824 

-10.362 

-  7.89 

.319 

C  -  5 
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